PADE APPROXIMANTS FOR FUNCTIONS WITH BRANCH POINTS - 
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Abstract. Let / be a germ of an analytic function at infinity that can be analytically continued 
along any path in the complex plane deprived of a finite set of points, / £ A(C \ A), J) A < 00 • 
J. Nuttall has put forward the important relation between the maximal domain of / where 
the function has a single-valued branch and the domain of convergence of the diagonal Pade 
approximants for /. The Pade approximants, which are rational functions and thus single- 
valued, approximate a holomorphic branch of / in the domain of their convergence. At the same 
time most of their poles tend to the boundary of the domain of convergence and the support 
of their limiting distribution models the system of cuts that makes the function / single-valued. 
Nuttall has conjectured (and proved for many important special cases) that this system of cuts 
has minimal logarithmic capacity among all other systems converting the function / to a single- 
valued branch. Thus the domain of convergence corresponds to the maximal (in the sense of 
minimal boundary) domain of single- valued holomorphy for the analytic function / £ _4(C\A). 
The complete proof of Nuttall's conjecture (even in a more general setting where the set A has 
logarithmic capacity 0) was obtained by H. Stahl. In this work, we derive strong asymptotics for 
the denominators of the diagonal Pade approximants for this problem in a rather general setting. 
We assume that A is a finite set of branch points of / which have the algebro-logarithmic character 
and which are placed in a generic position. The last restriction means that we exclude from our 
consideration some degenerated "constellations" of the branch points. 



1. Introduction 

Let / be a function holomorphic at infinity. Then / can be represented as a power series 

OO n 

(i-i) = 

k=0 

A diagonal Pade approximant to / is a rational function [n/n]f = p n /q n of type (n,n) (i.e., 
deg(p„), deg(<7„) < n) that has maximal order of contact with / at infinity [42, 8]. It is obtained 
from the solutions of the linear system 

(1.2) Rn{z) := q n {z)f(z) - Pn (z) = O (l/z n+1 ) as z -> 00 

whose coefficients are the moments fk in (1.1). System (1.2) is always solvable and no solution of 
it can be such that q n = (we may thus assume that q n is monic). In general, a solution is not 
unique, but yields exactly the same rational function [n/n]f. Thus, each solution of (1.2) is of the 
form (lp n ,lq n ), where (p n ,qn) is the unique solution of minimal degree. Hereafter, (p ra ,g n ) will 
always stand for this unique pair of polynomials. 

Pade approximant [n/n]f as well as the index n are called normal if deg(<7„) = n [34, Sec. 2.3]. 
The occurrence of non-normal indices is a consequence of overinterpolation. That is, if n is normal 
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f(z) - [n/n] f (z) ~ z-( 2n+/+1 ) as z -> oo 



index and 1 

' J/l 

for some I > 0, then [n/n]f = [n + j/n + j]f for j € {0, . . . , £}, and + Z + 1 is normal. 

Assume now that the germ (1.1) is analytically continuable along any path in C \ A for some 
fixed set A. Suppose further that this continuation is multi- valued in C\^4, i.e., / has branch-type 
singularities at some points in A. For brevity, we denote this by 

(1.3) feA(C\A). 

The theory of Pade approximants to functions with branch points has been initiated by J. Nuttall. 
In the pioneering paper [35] he considered a class of functions (1.3) with an even number of branch 
points (forming the set ^4) and principal singularities of the square root type. Convergence in 
logarithmic capacity [45, 46] of Pade approximants, i.e., 

(1.4) Ve > 0, lim cp ({z e K : \f(z) - [n/n] f {z)\ > e}) = 0, 

71— >OC 

was proven uniformly on compact subsets of C\ A, where A is a system of arcs which is completely 
determined by the location of the branch points. Nuttall characterized this system of arcs as a 
system that has minimal logarithmic capacity among all other systems of cuts making the function 
/ single-valued in their complement. That is, 

(1.5) cp(A) = min cp(&D) , 

dD-.DEVf 

where we denoted by T>f the collection of all connected domains containing the point at infinity in 
which / is holomorphic and single- valued. 

In that paper he has conjectured that for any function f in A(C \ A) with any finite number of 
branch points that are arbitrarily positioned in the complex plane, i.e., 

(1.6) jjA<oo and AcC, 

and with an arbitrary type of branching singularities at those points, the diagonal Pade approxi- 
mants converge to f in logarithmic capacity away from the system of cuts A characterized by the 
property of minimal logarithmic capacity. 

Thus, Nuttall in his conjecture has put forward the important relation between the maximal 
domain where the multi- valued function / has single- valued branch and the domain of convergence 
of the diagonal Pade approximants to / constructed solely based on the series representation (1.1). 
The Pade approximants, which are rational functions and thus single- valued, approximate a single- 
valued holomorphic branch of / in the domain of their convergence. At the same time most of 
their poles tend to the boundary of the domain of convergence and the support of their limiting 
distribution models the system of cuts that makes the function / single- valued (see also [36]). 

The complete proof of Nuttall's conjecture (even in a more general setting) was taken up by 
H. Stahl. In a series of fundamental papers [47, 48, 49, 50, 52] for a multi-valued function / e 
A(C \ A) with cp(.A) = (no more restrictions!) he proved: the existence of a domain D* € T>f 
such that the boundary A = satisfies (1.5); weak (n-th root) asymptotics for the denominators 
of the Pade approximants (1.2) 

(1.7) lim - log \q n (z)\ = -V u *(z), z e D* , 

n— foo ri 

where V UA := — / log \z — t\ duj&(t) is the logarithmic potential of the equilibrium measure wa, 
minimizing the energy functional /(/x) := J V fl (z) dfi(z) among all probability measures /i on A, 
i.e., I{u>a) '■= m i n ^(A)=i ^(m); convergence theorem (1.4). 

■'"We say that u(z) ~ b(z) if < lim inf z _ ! . 0o \a(z)/b(z)\ < limsup z _ >00 \a(z)/b(z)\ < oo. 
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The aim of the present paper is to established the strong (or Szego type, see [55, Ch. XII]) 
asymptotics of the Nuttall-Stahl polynomials q n . In other words, to identify the limit 

lim |L =? in D*, 

where the polynomials q n are the denominators of the diagonal Pade approximants (1.2) to func- 
tions (1.3) satisfying (1.6) and $ is a properly chosen normalizing function. 

Interest in the strong asymptotics comes, for example, from the problem of uniform convergence 
of the diagonal Pade approximants. Indeed, the weak type of convergence such as the convergence 
in capacity in Nuttall's conjecture and Stahl theorem is not a mere technical shortcoming. Indeed, 
even though most of the poles (full measure) of the approximants approach the system of the 
extremal cuts A, a small number of them (measure zero) may cluster away from A and impede 
the uniform convergence. Such poles are called spurious or wandering. Clearly, controlling these 
poles is the key for understanding the uniform convergence. 

There are many special cases of the Nuttall-Stahl polynomials that have been studied in detail 
including their strong asymptotics. Perhaps the most famous examples are the Pade approximants 
to functions 1/v ' z 2 — 1 and \ z 2, — 1 — z (the simplest meromorphic functions on a two sheeted 
Riemann surface of genus zero) where the Nuttall-Stahl polynomials q n turn out to be the classical 
Chebyshev polynomials of the first and second kind, respectively. The study of the diagonal Pade 
approximants for functions meromorphic on certain Riemann surfaces of genus one by means of 
elliptic functions was initiated in the works of S. Duma [19] and N.I. Akhiezer [2], see also [33] by 
E.M. Nikishin. Supporting his conjecture, Nuttall considered two important classes of functions 
with branch points for which he obtained strong asymptotics of the diagonal Pade approximants. 
In a joint paper with S.R. Singh [41], a generalization of the class of functions considered in [35] 
(even number of quadratic type branch points) was studied. Peculiarity of this class as well as 
its prototype from [35] is that A, the system of extremal cuts (1.5), consists of non-intersecting 
analytic arcs, see Figure 1A. 





(A) (B) 

FIGURE 1. The left-hand figure depicts the case where A is comprised of non- 
intersecting analytic arcs, and the right-hand figure illustrates the case where three arcs 
share a common endpoint. 



In the paper [39] (see also [37]) Nuttall investigated the behavior of the Pade approximants for 
functions with three non-collinear branch points. Namely, functions of the form 
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(1.8) f{z) =H(z e.,r , <*i e C : 5>,- = . 



=i 



Analytic arcs of the system of extremal cuts for these functions (contrary to the functions from 
the previous class) share a common endpoint, see Figure IB. In order to shed some light on the 
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behavior of the spurious poles, Stahl studied strong asymptotics of the diagonal Pade approximants 
for hyperelliptic functions [51], see also [52]. 

An important feature of the diagonal Pade approximants, which plays a key role in the study 
of their asymptotics, is the orthogonality of their denominators. It is quite simple to see that (1.2) 
and the Cauchy theorem (taking into account the definition of T>f in (1.5)) lead to 



q n (z) z 3 f(z) dz = o, j e {o, ... ,n- 1}, 

dD:DeV s 

where the integral is taken along the orientated boundary dD of a domain D from T>f. For the 
extremal domain D* £ T)f of / satisfying (1.3) and (1.6), the boundary A = dD* consists of a 
finite union of analytic Jordan arcs. Hence, choosing an orientation of A (as a set of Jordan arcs), 
we can introduce in general complex-valued weight function 

(1.9) P (t) = (f+-f-)(t), teA, 

which turns q n into non-Hermitian orthogonal polynomials. That is, 

(1.10) J q n (t)t j p(t)dt = 0, j G {0, ... ,71 — 1}. 

A 

Asymptotic analysis of the non-Hermitian orthogonal polynomials is a difficult problem substan- 
tially different from the study of the asymptotics of the polynomials orthogonal with respect to a 
Hermitian inner product, i.e., the case where p is real- valued and Acl. 

In [49], Stahl developed a new method of study of the weak (n-th root) asymptotics (1.7) of the 
polynomials orthogonal with respect to complex-valued weights. As discussed above, this resulted 
in the proof of Nuttall's conjecture. The method of Stahl was later extended by A. A. Gonchar 
and E.A. Rakhmanov in [28] to include the weak (n-th root) asymptotics of the polynomials 
orthogonal with respect to varying complex- valued weights, i.e., to include the case where the weight 
function p := p n in (1.10) depends on n, the degree of the polynomial q n . Orthogonal polynomials 
with varying weights play an important role in analysis of multipoint Pade approximants, best 
(Chebyshev) rational approximants, see [27, 26], and in many other applications (for example in 
description of the eigenvalue distribution of random matrices [14]). 

The methods of obtaining strong asymptotics of the polynomials orthogonal with respect to a 
complex weight are based on a certain boundary-value problem for analytic functions (Riemann- 
Hilbert problem). Namely, 

(1.11) R+-R-=q nP on A, 

where R n , defined in (1.2), are the reminder functions for Pade approximants (or functions of the 
second kind for polynomials (1.10)), which also can be expressed as 

(1.12) m-jjmp.*., Jtw-o^) - »-.». 

The boundary-value problem (1.11) naturally follows from (1.2) and the Sokhotskii-Plemelj for- 
mulae. This approach appeared in the works of Nuttall in connection with the study of the strong 
asymptotics of the Hermite-Pade polynomials, see the review [38] . In [40] Nuttall transformed the 
boundary condition (1.11) into a singular integral equation and on this basis obtained the formulae 
of strong asymptotics for polynomials (1.10) orthogonal on the interval A := [—1, 1] with respect 
to a holomorphic complex-valued weight 

p(x) peH(A), p^Oon A:=[-l,l], 

VI — x z 
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where if (A) is a class of functions holomorphic in some neighborhood of A. Here p can also be a 
complex- valued non- vanishing Dini-continuous function on [—1,1] [11]. The most general known 
extension of this class of orthogonal polynomials is due to S.P. Suetin [53, 54] who considered 
the convergence domain D* £ T>f for the function 1/ *J (t — e±) ■ ■ ■ (t — e2 g +2), &j € C, when the 
boundary A = dD* consists of g + 1 disjoint Jordan arcs (like in [41], see Figure 1A). Elaborating 
on the singular integral method of Nuttall, he derived strong asymptotics for polynomials (1-10) 
orthogonal on A with respect to the complex weight 



^/{t - ei) • • • (t - e 2 g+2) 

where p is a Holder continuous and non- vanishing function on A. In [10], L. Baratchart and the 
second author have studied strong asymptotics for polynomials (1-10) via the singular integral 
method in the elliptic case g = 1, but under the assumptions that p is Dini-continuous and 
non- vanishing on A, while the latter is connected and consists of three arcs that meet at one of 
the branch points (exactly the same set up as in [39], see Figure IB). The strong asymptotics 
of Nuttall-Stahl polynomials arising from the function (1.8) was derived in the recent work [21] 
in three different ways, including singular integral equation method of Nuttall and the matrix 
Ricmann-Hilbert method. 

The latter approach facilitated substantial progress in proving new results for the strong asymp- 
totics of orthogonal polynomials and is based on a matrix-valued Riemann-Hilbcrt boundary value 
problem. The core of the method lies in formulating a Ricmann-Hilbert problem for 2x2 matri- 
ces (due to Fokas, Its, and Kitaev [22, 23]) whose entries are orthogonal polynomials (1-10) and 
functions of the second kind (1-12) to which the steepest descent analysis (due to Deift and Zhou 
[17]) is applied as n —¥ oo. This method was initially designed to study the asymptotics of the 
integrable PDEs and was later applied to prove asymptotic results for polynomials orthogonal on 
the real axis with respect to real-valued analytic weights, including varying weights (depending 
on n) [16, 15, 31, 32] and related questions from random matrix theory. It also has been noticed 
[7, 3, 29, 12] (see also recent paper [13]) that the method works for the non-Hermitian orthogonality 
in the complex plane with respect to complex-valued weights. 

In the present paper we apply the matrix Riemann-Hilbert method to obtain strong asymptotics 
of Pade approximants for functions with branch points (i.e., we obtain strong asymptotics of 
Nuttall-Stahl polynomials). To capture the geometry of multi-connected domains we use the 
Riemann theta functions as it was done in [16], but keep our presentation in the spirit of [4, 6]. 

This paper is structured as follows. In the next section we introduce necessary notation and 
state our main result. In Sections 3 and 4 we describe in greater detail the geometry of the problem. 
Namely, Section 3 is devoted to the existence and properties of the extremal domain D* for the 
functions of the form (1.3) and (1.6). Here, for completeness of the presentation, we present some 
results and their proofs from the unpublished manuscript [43]. Section 4 is designed to highlight 
main properties of the Riemann surface of the derivative of the complex Green function of the 
extremal domain D* . Sections 5 and 6 are devoted to a solution of a certain boundary value 
problem on 91 and are auxiliary to our main results. In the last three sections we carry out the 
matrix Riemann-Hilbcrt analysis. In Section 7 we state the corresponding matrix Riemann-Hilbcrt 
problem, renormalize it, and perform some identical transformations that simplify the forthcoming 
analysis. In Section 8 we deduce the asymptotic (as n — » oo) solution of the initial Ricmann-Hilbert 
problem and finally in Section 9 we derive the strong asymptotics of Nuttall-Stahl polynomials. 
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2. Main Results 



The main objective of this work is to describe the asymptotics of the diagonal Pade approximants 
to algebraic functions. We restrict our attention to those functions that have finitely many branch 
points, all of which are of an integrable order, with no poles and whose contour of minimal capacity 
satisfies some generic conditions, Section 2.1. Such functions can be written as Cauchy integrals 
of their jumps across the corresponding minimal capacity contours and therefore we enlarge the 
considered class of functions to Cauchy integrals of densities that behave like the non-vanishing 
jumps of algebraic functions, Section 2.2. It turns out that the asymptotics of Pade approximants 
is described by solutions of a specific boundary value problem on a Riemann surface corresponding 
to the minimal capacity contour. This surface and its connection to the contour are described in 
Section 2.3, while the boundary value problem as well as its solution are stated in Section 2.4. The 
main results of this paper are presented in Section 2.5. 

2.1. Functions with Branch Points. Let / be a function holomorphic at infinity that extends 
analytically, but in a multi- valued fashion, along any path in the extended complex plan that omits 
finite number of points. That is, 



Without loss of generality we may assume that /(oo) = since subtracting a constant from / 
changes the Pade approximant in a trivial manner. 

We impose two general restrictions on the functions (2.1). The first restriction is related to 
the character of singularities at the branch points. Namely, we assume that the branch points are 
algebro-logarithmic. It means that in a small enough neighborhood of each a k the function / has 
a representation 



where —1 < a(a,k) < and h\,h2 are holomorphic around a k - The second restriction is related 
to the disposition of the branch points. Denote by D* the extremal domain for / in the sense of 
Stahl. It is known, Proposition 8, that 



where 1J is a finite union of open analytic Jordan arcs and E is a finite set of points such that 
each element of E is an endpoint for at least one arc A^, Figure 2. In what follows, we suppose 



(2.1) 



/ G A(C \ A) , 



A := {a k }, 2 < flA < oo. 



(2.2) 




(2.3) 





FIGURE 2. Schematic example of A, depicting the set E — {a\, . . . , a^} U {6i, 62}, the 
arcs A*,, and their orientation. 



that the points forming A are in a Generic Position (GP). 

Condition GP. We assume that 

(i) each point in E P| A is incident with exactly one arc from [J A^ ; 
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(ii) each point in E \ A is incident with exactly three arcs from (J A&. 

The above condition describes a generic case for the set A. Meaning that if the set A does not 
satisfy this condition, then there is a small perturbation of the position of these points such that 
new set A obeys Condition GP. 

Denote by <?a the Green junction for D* with a pole at infinity, Section 3. That is, g A is the 
unique function harmonic in D* \ {00} having zero boundary values on A that diverges to infinity 
like log \z\ as \z\ — > 00. It is known [45, Thm. 5.2.1] that the logarithm of the logarithmic capacity 
of A is equal to 

logcp(A)= lim (\og\z\- g A (z)). 
As shown in [43], see also (3.9) further below, it holds that 

(2.4) h(z) := (2d z9A )(z) = - + ••• = 

z 

where 2d z := d x — id y , B is a monic polynomial of degree deg(A) — 2, and 

p 

A(z) := Y[(z-a k ), {a u ...,a p } :=AnE. 

k=l 

Since g A = on A, so is its tangential derivative at each smooth point of A + U A - . Hence, h(t)rt, 
t G A*, is purely imaginary, where r t is the complex number corresponding to the tangent vector 
at t to A. In particular, the integral of h(t)dt along is purely imaginary. 

Condition GP has the following implications on B [44, Section 8] . Let m be the number of the 
connected components of A. Then B has p — 2m simple zeros that we denote by b%, . . . , 6 p _2m and 
all the other zeros are of even multiplicities. In particular, 

(2.5) E := {oi, . . . ,a p } U {b%, . . . ,b p ^ 2m }- 
If we set g := p — m — 1, then \E\ = 2g + 2. Moreover, we can write 

p— 2m g 

(2.6) b{z)= n (* - h i) n ( z -^) 2 

j — l j—p—2m+l 

where the elements of {6 p -2m+ij • ■ • , b g } are the zeros of B of even multiplicities , each listed 
m(,./2 times. 

2.2. Cauchy-type Integrals. Let / be a function of the form (2.1)-(2.2) with contour A in 
(2.3) satisfying Condition GP. We orient the arcs A^ comprising A so that the arcs sharing a 
common endpoint either are all oriented towards this endpoint or away from it, see Figure 2. As 
the complement of A is connected, i.e., A forms no loop, such an orientation is always feasible and, 
in fact, there are only two such choices which are inverse to each other. According to the chosen 
orientation we distinguish the left (+) and the right (— ) sides of each arc. Then 

P-7) fw-fSt^n®* .en; 



j a t — z 2ni ' 
where the integration on A is taking place according to the chosen orientation. 

Let e G E \ A. Then e is incident with exactly three arcs, which we denote for convenience 
by A ej -, j G {1,2,3}. Since e is not a point of branching for /, the jumps (/ + — /~)|a c k are 
holomorphic around e and enjoy the property 

(2.8) (/+ - /-)| Ae , 1 + (/+ - r K, a + (/ + - .r)iA B , 3 = o, 

where we also used the fact that the arcs A e j have similar orientation as viewed from e. 
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Set a(e) := for each e E E \ A and fix a, b € E that are adjacent to each other by an arc 
Afe C A. Then the jump of / across A^ can be written as 

(2.9) (f+-f-)(z)=w k (f;z)(z-a) a ^(z-b) a ^, zeA fe , 

where we fix branches of (z — a) a ^ and (z — b) a ^ that are holomorphic across A/, and w k (f] ■) 
is a holomorphic and non- vanishing function in some neighborhood of A&. 

Keeping in mind (2.8) and (2.9), we introduce the following class of weights on A. 

Definition 1. A weight function p on A belongs to the class Wa if 

(2.10) jO| Ae>1 + P|A e>2 + P|A e ,3 = 

in a neighborhood of each e € E \ A, where A e j 7 j E {1, 2, 3}, are the arcs incident with e; and 

(2.11) p, Afc (z) - w k {z){z-a) a °{z-b) ab 

on A k , incident with a,b, where w k is holomorphic and non-vanishing in some neighborhood of 
Afc and {(z — e) a '} e£ E is a collection of functions holomorphic in some neighborhood of A \ {e}, 
a e > — 1 and a e = for e E E\ A. 

For a weight p E Wa, we set 

(2.12) m:=f z 6 D* ' 

JA t — Z ITTl 

It follows from (2.7)-(2.9) that a function / satisfying (2.1)-(2.2) and Condition GP, and whose 
jump /+ — / _ is non-vanishing on A \ {a\, . . . , a p } can be written in the form (2.12) for a Wa- 
weight. 

2.3. Riemann Surface. Denote by £R the Riemann surface of h defined in (2.4). We represent 
£K as a two-sheeted ramified cover of C constructed in the following manner. Two copies of C are 
cut along every arc A k . These copies are joined at each point of E and along the cuts in such a 
manner that the right (resp. left) side of the arc A k belonging to the first copy, say 0v-°>, is joined 
with the left (resp. right) side of the same arc A k only belonging to the second copy, DV* 1 '. It can 
be readily verified that 9\ is a hyperelliptic Riemann surface of genus g. 

According to our construction, each arc A k together with its endpoints corresponds to a cycle, 
say Lfc, on 9L We set L := [j k L k , denote by ir the canonical projection ir : £H — > C, and define 
D (k) ._ m (k) n yr-ip*) and z lk) ._ D (k) n 

for fc E {0, 1} and z g D* . We orient each L k in such a manner that remains on the left when 
the cycle is traversed in the positive direction. For future use, we set }^ =1 to be such point 
on *H that 

(2.13) Tr{bf ] )=b 3 and bfeD^UL, je{l,...,g}. 

For any (sectionally) meromorphic function r on £H we keep denoting by r the pull-back function 
from onto C and we denote by r* the pull-back function from D^ 1 ) onto C. We also consider 
any function on C naturally defined on JH'°). In particular, ft, is a rational function over 1H such 
that h* = —h (as usual, a function is rational over 9\ if the only singularities of this function on 9i 
are polar). 

Denote by {a fc }^ =1 and {bk} k =i tne following homology basis for d\. Let Cj(A), j G {1, . . . , m}, 
be the connected components of A. Set pj := \{ax, ■ ■ ■ , a p }(lCj(A)\. Clearly, p = Y^7=\Pr Relabel, 

if necessary, the points {ai, . . . , a p } in such a manner that {a Po _| ^ Pj _ 1 +i, ■ • ■ , a po ^ ^ Pj } C Cj(A), 

where p Q := 0. Then for each j E {2, . . . , m} we choose pj — 1 b-cycles as those cycles Lf. that 
contain the points a poH hPj-1+2, ■ • • , a p A Vvi > anc ^ we c h° ose Pi~ 2 b-cycle as those cycles L k that 
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contain 03, ... , a Pl . We assume that the orientation of the b-cycles is induced by the orientation of 




FIGURE 3. The choice of the b-cycles and the parts of the a-cycles belonging to 
(thicker lines necessarily oriented towards the corresponding b-cycles). 



the corresponding cycles L k . The a-cycles are chosen to be mutually disjoint except at et2, which 
belongs to all of them. It is assumed that each cycle a k intersect the corresponding cycle b fc only 
at one point, the one that belongs to {03, . . . , a p }, and that 

A a k := n(a k n D (0) ) = n{a k n £> (1) ). 

The a-cycles are orientated in such a manner that the tangent vectors to a^b^ form the right 
pair at the point of their intersection. We also assume that each arc A^ naturally inherits the 
orientation of afe n D^°\ In particular, the + side of a^ n and the — side of a k H project 
onto the + side of A?, see Figure 3. We set 



m : 



m \ (J (afc U b fe ) and m := <K \ Q a k 



k=l 



k=l 



(observe that 91 is a simply connected subdomain of 5H). 
Define (see also Section 4.2) 



(2.14) 



$(z) 



exp 



h(t)dt 
01 J 



for z € m. 



Then $ is a holomorphic and non- vanishing function on 91 except for a simple pole at 00 ^ and a 
simple zero at 00 W whose pull-back functions are reciprocals of each other, i.e., 



(2.15) 



1 in D a :=D*\{jAl 



k=l 



Furthermore, $ possesses continuous traces on both sides of each a- and b-cycle that satisfy 

$+ _ J exp{27T^ fc } on a k , 
$- \ exp{27rir fc } on b fe , 

where the constants io k and r k are real and can be expressed as 
(2.17) u k 



■ — <f> h(t)dt and r k := — !— (p h(t)dt, 
2™ / b , 2m J 
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k € {1, . .. ,g}. In fact, it holds that uj k = uj&('K(Lk)), where u>a is the equilibrium measure of A 
[45]. Moreover, it is true that 

(2.18) $ (z) = _i_ + ( i) as z^oo, |£| = 1. 

In what follows, we shall assume without loss of generality that £ = 1. Indeed, if £ ^ 1, set 
:= {£z : z 6 A} and p^(z) := p(£,z), z € Ag, where p is a function defined in (2.11). Then 

Pi( z ) K&) and i n / n ]h( z ) = \- n / n b(£, z )' 2GC\A ? . 
Thus, the asymptotic behavior of [n/n]^ is entirely determined by the asymptotic behavior of 
[n/n]p v Moreover, it holds that $a 5 (z) = $a(£z) an d therefore $a 5 (z) = z/cp(A) + 0(1). That 
is, we always can rotate the initial set up of the problem so that (2.18) holds with £ = 1 without 
altering the asymptotic behavior. 

Recall that a Riemann surface of genus g has exactly g linearly independent holomorphic dif- 
ferentials (see Section 4.1). We denote by 

dfi := (dfl 1 ,...,dQ g ) T 
the column vector of g linearly independent holomorphic differentials normalized so that 

(2.19) W dQ = e k for each k e {1, . . . , g}, 

J afc 

where {e k } 9 k=1 is the standard basis for M. 9 and (F is the transpose of e. Further, we set 

(2.20) Bn := <f dfl k 

It is known that Bq is symmetric and has positive definite imaginary part. 



2.4. Auxiliary Boundary Value Problem. Let p e Wa- Define 
(2.21) 



u := (uj 1 ,...,u! g ) T , 
t := (ri,...,T ff ) , 



for some fixed determination of log(p/h + ) continuous on A \ E, where the constants ujj and Tj 
were defined in (2.17) and we understand that \og(p/h + ) on L is the lift \og(p/h + ) o n. 

Further, let {tj} be an arbitrary finite collection of points on SH. An integral divisor corre- 
sponding to this collection is defined as a formal symbol ^ tj . We call a divisor ^2 tj special if it 
contains at least one pair of involution-symmetric points; that is, if there exist tj =/= t k such that 
7r(tj) = 7r(tfc) or multiple copies of points from E (with a slight abuse of notation, we keep using 
E for ^(E)). 

Given constants (2.21) and points (2.13), there exist divisors J2j=i t n .j, see Sections 4.3 and 6.1 
further below, such that 

(2.22) / dfi = c p + n(Cu + But) (mod periods dfi] , 

where the path of integration belongs (for definiteness, we shall consider each endpoint of 
integration belonging to the boundary of IK as a point on the positive side of the corresponding a- 
or b-cycle) and the equivalence of two vectors c , e e C 9 is defined by c = e ^mod periods dtfj if 
and only if c — e = j + B^fh for some j, fa £ Z 9 . 
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Proposition 1. Solutions of (2.22) are either unique or special. If (2.22) is not uniquely solvable 
for some index n, then all the solutions for this index assume the form 

g-2k k 

£t,. + £(z<°> + 

2=1 2=1 

where the divisor X)f=i fc * s fixed and non-special and {zj}j_-y are arbitrary points in C. 
If for some index n a divisor solving (2.22) has the form 

g-i 

5^ti + fcoo<°> + (I - k)oo^ 

i=l 

with I > 0, k £ {0, . . . , I}, and non-special X)f=i such that |7r(tj)| < oo, then 

g-i 

+ j)oo {0) + (l-k- j)oo (1) 

i=l 

solves (2.22) for the index n + j for each j £ {~k, . . . ,1 — k}. In particular, (2.22) is uniquely 
solvable for the indices n — k and n + I — k. 

If y^j—i t n j uniquely solves (2.22) and does not contain oo^ k \ k £ {0, 1}, then (2.22) is uniquely 
solvable for the index n — (— l) fc and {t„.j}| =1 f~l {t n -(-i) fc j = ■ 

Remark 1.1. Propositions 1 says that the non-unique solutions of (2.22) occur in blocks. The last 
unique solution before such a block consists of a non-special finite divisor and multiple copies of 
oo^ . Trading one point oo^ 1 ) for oo^ ) and leaving the rest of the points unchanged produces a 
solution of (2.22) (necessarily non-unique as it contains an involution-symmetric pair oo' 1 ) -t-oo' -') 
for the subsequent index. Proceeding in this manner, a solution with the same non-special finite 
divisor and all the remaining points being oo(°) is produced, which starts a block of unique solutions. 
In particular, there cannot be more than g ~ 1 non-unique solutions in a row. 

Definition 2. In what follows, we always understand under 53f=i ^n,j either the unique solution 
of (2.22) or the solution where all the involution-symmetric pairs are taken to be oo' 1 ) + oo^ -*. 
Under this convention, given e > 0, we say that an index n belongs to N £ C N if and only if 

(i) the divisor 53f=i ^"n,j satisfies |*7r(t n j)| < 1/e for all t„j £ £HW; 

(ii) the divisor X}f=i*n-ij satisfies |7r(t n _ij)| < 1/e for all t n -\ t j £ . 

To show that Definition 2 is meaningful we need to discuss limit points of {J2j=i tn,j}neW, 
N' C N, where convergence is understood in the topology of 9\ 9 /T, g , 9l 9 quotient by the symmetric 
group £ 9 . The following proposition shows that these limiting divisors posses the same block 
structure as the divisors themselves. 

Proposition 2. Let N' be such that all the limit points of { 2f=i *n.i} TieN / assume the form 

g—2k—lo—li k 

(2.23) £ t, + £ (zP + z™) + W 0) + hoc*™ 

i=l i=l 

for a fixed non-special divisor 53f=i ° 1 < °°j an d arbitrary {zi}^ =1 C C. Then all 

the limit points of the sequence { X)f=i ^n+j-i} n ew > 3 ^ {~^o — k, . . . ,h + k}, assume the form 

g—2k—lo—h k' 

(2.24) £ *< + E H 0) + ™i X) ) + («o + J + * - fc')°° (0) + Oi - J + fc - 
i=l i=l 
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where < k! < min {Iq + k + j, h + fc — j} and {wi}k =1 C C. 

// { Ylj—i tn,i} neN , converges to a non-special divisor 5Zf=i that does not contain oo' fc \ k € 
{0,1}, £/ien £/ie sequence { 5Zi=i tn-(-i) fc ,i} rlGN / a ' so converges, say to X)f=i w j'j which is non- 
special, and {tj}j =1 H {wj}| =1 = 0. 

Remark 2.1. Proposition 2 shows that the sets N £ arc well-defined for all e small enough. Indeed, 
let { J2i=i t n ,i} neW be a subsequence that converges in d\ 9 /S s (it exists by compactness of 9^). 
Naturally, the limiting divisor can be written in the form (2.23). Then it follows (2.24) that the 
sequence { Yn=i *n-h-k,i} ngN , converges to YfiZi ~ °~ 1 *» + (h + lo + 2fc) oo^. Further, by the 
second part of the proposition, the sequence { 5Zf=i ^n-l -k-i,i} ne ^, also converges and the limit, 
sa y zCf=i w i: does not contain oc/ 1 ). Thus, {n — la — k : n € N'} C N £ for any e satisfying 
e|7r(ti)| < 1 if U e <H (0) and e|7r(w. i )| < 1 if w, <= 

Equipped with the solutions of (2.22), we can construct the Szego functions of p on 9t, which 
are the solutions of a sequence of boundary value problems on L. 

Proposition 3. For each n € N there exists a function, say S n , with continuous traces on both 
sides of (L U (Jfe=i a fe) \ E such that S n & n is meromorphic in Vl \ L and 

(2.25) (5 n $™)" = (p/h+)(S n $ n )+ on L\E. 

If we let m(t) to be the number of times, possibly zero, t appears in {t n ,j}| =1 , then S n is non- 
vanishing and finite except for 

( \S n (zW)\ ~\Z- a |™(a)/2-(-l) fc (l+2a )/4 flg z (k) _^ ffl g { aj } P j=v 

(2.26) J |S„(z«)| ~ | z _ h |m(6)/2-l/2+(-l)V4 as z (fc) {6 J -}JlJ m , 

( \S n {zM)\~\z-bF«>W as ^^be{bf^ =p _ 2m+v 

and has a zero of multiplicity m(t) at each t £ {t n .j}j =1 \ ( { a j}^ = i U an d m b/2 is the 

multiplicity of b in {bj} 9 j=p _ 2m+1 - 

Conversely, if for given neN t/iere exists a function S with continuous traces on (L U Ufe=i a fc)\ 
E such that S3> n is meromorphic in and S satisfies (2.25) and (2.26) with Y]j—i ^n.j replaced 

by some divisor J]tj, then Y]tj solves (2.22) /or the index n and S = pS n for a polynomial p 
such that (p on) — ^ tj — ^n,j ■ 

Finally, given e > and e > 0, i/iere exists constant C e . £ < oo such that 

S„_i(z) S»(oo(°>) 



(2.27) 



S n («) S„-i(oc«) 

/or n G N £ and z 6 £H n , c := 9t\ U" =1 A^(t„j), where N e (t n j) is a connected neighborhood oft n j 
such that iv(N e (t n j)) is the e-ball centered at t n j in the spherical metric 2 . 

Remark 3.1. The integers m(e), e e E, in the first two lines of (2.26) are either or 1 as otherwise 
J2j=i tnj would be special. 

Remark 3.2. The estimate in (2.27) cannot be improved in a sense that if for some subsequence of 
indices e n — > 0, l/e n := max { max tn _ 1 ^g^ 1 ) (""(tn-ij)! , m ax tn . e s>t(o) |7r(t ni j)| }, then C £i£ll — > oo. 

Remark 3.3. We would like to stress that S n is unique for n € N e , £ > 0, as (2.22) is uniquely 
solvable for all such indices. 

Propositions 1-3 are proved in Section 6. 

2 That is, dist(z 1 ,« 2 ) = 2|*i - ^ 2 |(1 + |^i| 2 )- 1 / 2 (l + |^ 2 | 2 )~ 1/2 if |*ll.l«2| < °o and dist(z,oo) = 2(1 + \z\ 2 )- 1 / 2 . 
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2.5. Main Theorem. Let {[n/n]p} n ^ be the sequence of diagonal Pade approximants to the 
function p. As before, denote by q n the denominator polynomial of [n/n]p (Nuttall-Stahl orthogonal 
polynomial (1.10)) and by R n the reminder function of [n/n]^ (1.2) (the function of the second 
kind (1.12) for q n ). Recall that by S n and 5* we denote the pull-back functions of S n on 91 from 
Z)(°) and to D*, respectively. Then the following theorem holds. 

Theorem 4. Let A be a minimal capacity contour as constructed in Section 2.1 subject to Condi- 
tion GP and assumption £ = 1 in (2.18). Further, let p be given by (2.12) with p € Wa, N e be as 
in Definition 2 for fixed e > 0, and S n be given by Proposition 3. Then for all n£N E it holds that 

In = (1 + V n i) J n S n <$> n + U n , 2 7*S' n _ 1 $ n - 1 , 



(2.28) 



R n = ( l + V nl)ln^ r +V n2ln -— r 



locally uniformly in D* , where \v n j\ < c{e)/n in C while v n j(oo) = and 

cp(A)" , . cp(A)"+ 1 
5„(oo) ^-1(00) 

Moreover, it holds locally uniformly in A \ E that 
(2.29) 

Qn = (1 + v nl ) 7 „ ((5„$") + + (£„* n )-) + W „27: ((5 n _ 1 $"" 1 ) + + ^-x^"- 1 )") , 



Before proceeding, we would like to make several remarks regarding the statement of Theorem 4. 

Remark 4.1. If the set A consists of two points, then A is an interval joining them. In this case 
the conclusion of Theorem 4 is contained in [38, 5, 32, 12]. Moreover, the Riemann surface 91 has 
genus zero and therefore <E> is simply the conformal map of D* onto {\z\ > 1} mapping infinity into 
infinity and having positive derivative there, while S n = S p is the classical Szego function. 

Remark 4.2. Notice that both pull-back functions S n and hS* are holomorphic D* . Moreover, S n 
has exactly g zeros on 91 that do depend on n. It can be deduced from (2.28) that q n has a zero in 
the vicinity of each zero of S n that belongs to D^°\ These zeros are called spurious or wandering 
as their location is determined by the geometry of 91 and, in general, they do not approach A with 
n while the rest of the zeros of q n do. On the other hand, those zeros of S n that belong to are 
the zeros of the pull-back function S* and therefore describe locations of the zeros of R n (points 
of overinterpolation) . 

Remark 4.3. Even though our analysis allows us to treat only normal indices that are also asymp- 
totically normal, formulae (2.28) illuminate what happens in the degenerate cases. If for an index 
n the solution of (2.22) is unique and contains I copies of oo'- 1 '', the function vanishes at infinity 
with order /. The latter combined with the second line of (2.28) shows that [n/n]p is geometrically 
close to ovcrintcrpolating p at infinity with order I. Then it is feasible that there exists a small 
perturbation of p (which leaves the vector c p unaltered) that turns the index n into a last normal 
index before a block of size I of non-normal indices, which corresponds to the fact that solutions 
of (2.22) are special for the next I — 1 indices and the solution for the index n + 1 contains I copies 
ofooW. 

Observe that p — [n/n]p — R n /q n by (1.2) applied with / := p. Thus, the following result on 
uniform convergence is a consequence of Theorem 4. 
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Corollary 5. Under the conditions of Theorem 4, it holds for neN E that 
(2.30) ? _[ n/n]?= [i + (i/ n )]g^. 

in D* H ^(Dt^e), where 0{l/n) is uniform for each fixed e > 0. 

3. Extremal Domains 

In this section we discuss existence and properties of the extremal domain D* G T>f for the 
function /, holomorphic at infinity that can be continued as a multi- valued function to the whole 
complex plane deprived of a polar set A, see (1.3). Recall that the compact set A := 3D* defined 
in (1.5) makes / single-valued in its complement and has minimal logarithmic capacity among all 
such compacta. 

As mentioned in the introduction, the question of existence and characteristic properties of A 
was settled by Stahl in the most general settings. Namely, he showed that the following theorem 
holds [47, Theorems 1 and 2] and [48, Theorem 1]. 

Theorem S. Let f e A(C\A) with cp(A) = 0. Then there exists unique D* e V f , A = C\D*, 
the extremal domain for f , such that 

cp(A) < cp(frD) for any DeV f , 

and if cp(A) = cp(dD) for some D e Vf, then D C D* and cp(D* \ D) =0. Moreover, 
A := Eq U Ei U 1J Afc, where Eq C A, Ei is a finite set of points, and A& are open analytic Jordan 
arcs. Furthermore, it holds that 

3.1) ^ — = ^ — on MA fc , 

where <?a is the Green function for D* and n. are the one-sided normals on each Aft. 

Let now / and A be as in (2.1). Denote by K, the collection of all compact sets K such that K 
is a union of a finite number of disjoint continua each of which contains at least two point from A 
and C\ K e V f . That is, 

{q<oo 
K : K= |J Kj, %(AnKj)> 2; Kj \ dKj = 0; K j DK l ^0, i^j; C \ K e V f 

Observe that the inclusion {C \ K : K G K,} C Df is proper. However, it can be shown using 
the monodromy theorem (see, for example, [9, Lemma 8]) that A E JC. Considering only functions 
with finitely many branch points and sets in JC allows significantly alter and simplify the proof of 
Theorem S, [43, Theorems 2 and 3]. Although [43] has never been published, generalizations of 
the method proposed there were used to prove extensions of Theorem S for classes of weighted 
capacities, see [30], [20] and [9]. Below, in a sequence of propositions, we state the simplified 
version of Theorem S and adduce its proof as devised in [43] solely for the completeness of the 
exposition. 

Proposition 6. There exists A € K. such that cp(A) < cp(K) for any K € K. 
Proof. Let {K n } be a sequence in JC such that 

cp(-fC„l — > inf cpf-ftT) =: c as n — > oo. 
Keic 

Then there exists R > such that K n C U)r := {z : \z\ < R} for all n large enough. Indeed, it is 
known [45, Theorem 5.3.2] that cp(A) > cp(7) > |diam(7), where 7 is any continuum in K n and 
diam(7) is the diameter of 7. As 7 contains at least two points from A, the claim follows. 



STRONG ASYMPTOTICS OF NUTTALL-STAHL POLYNOMIALS 



15 



For any K e JCr := JC H D fl and e > 0, set (K) t := {z : dist(z, K) < e}. We endow JCr with the 
Hausdorff metric, i.e., 

d H (K u K 2 ) := inf{e : K\ C (K 2 ) ei K 2 C (X x ) e }. 

By standard properties of the Hausdorff distance [18, Section 3.16], clos^^ {JCr), the closure of JCr 
in the dy-metric, is a compact metric space. Notice that a compact set which is the d#-limit 
of a sequence of continua is itself a continuum. Observe also that the process of taking the de- 
limit cannot increase the number of the connected components since the e-neighborhoods of the 
components of the limiting set will become disjoint as e — > 0. Thus, each element of clos^ (JCr) 
still consists of a finite number of continua each containing at least two points from A but possibly 
with multiply connected complement. However, the polynomial convex hull of such a set, that is, 
the union of the set with the bounded components of its complement, again belongs to JCr and 
has the same logarithmic capacity [45, Theorem 5.2.3]. 

Let A* € clos^ (JCr) be a limit point of {K n }. In other words, cLr (A*, K n ) — > as n — > oo, 
n e Ni C N. We shall show that 

(3.2) cp(A*) = c. 

To this end, denote by K e :— {z : gic(z) < log(l + e)}, e > 0, where gx is the Green function with 
pole at infinity for the complement of K. It can be easily shown [45, Theorem 5.2.1] that 

(3.3) cp(X £ ) = (1 + e)cp(X). 

Put Co := inf{cp(7)}, where the infimum is taken over all connected components 7 of K n and all 
n £ Nj. Recall that each component 7 of any K n contains at least two points from A. Thus, it 
holds that c > since cp(7) > |diam(7). 

We claim that for any e E (0, 1) and 5 < e 2 Co/2 we have that 

(3.4) (K n ) s c K 

for all n large enough. Granted the claim, it holds by (3.3) that 

(3.5) cp(A*) < (1 + e)cp(K„) 

since A* C (K n )g C K^. Thus, by taking the limit as n tends to infinity in (3.5), we get that 

(3.6) c< cp(A*) < (l + e)c, 

where the lower bound follows from the very definition of c since the polynomial convex hull of 
A*, say A, belongs to JC and has the same capacity as A*. As e was arbitrary, (3.6) yields (3.2) 
with A as above. 

It only remains to prove (3.4). We show first that for any continuum 7 with at least two points, 
it holds that 

e 2 

(3.7) dist( 7 ,7 e ) > -cp( 7 ), 

where 7 e := {z : g 1 (z) = log(l + e)}. Let f be a conformal map of {z : \z\ > 1} onto C \ 7, 
\&(oo) = 00. It can be readily verified that |5 , (z)z^ 1 | — > cp(7) as z — > 00 and that g 1 = log|<I'~ 1 |, 
where ^E 1-1 is the inverse of $ '. Then it follows from [25, Theorem IV. 2.1] that 

(3.8) |*'(3)|>cp( 7 )(i_ J_), \ z \>i. 

Let z\ £ 7 and z 2 £ 7 e be such that dist(7,7 e ) = \z\ — z 2 \. Denote by [zi, z 2 \ the segment joining 
Z\ and z 2 . Observe that & maps the annular domain bounded by 7 and j e onto the annulus 
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{z : 1 < \z\ < 1 + e}. Denote by S the intersection of ^ ((zi, £2)) with this annulus. Clearly, the 
angular projection of S onto the real line is equal to (1. 1 + e). Then 

dist( 7 ,7 e ) = / \dz\= [ \*'{z)\\dz\ > cp( 7 ) / (i_J_W| 



\ 1*1 / •'(1,1+e) \ 1*1 / J- 1" fc 

where we used (3.8). This proves (3.7) since it is assumed that e < 1. 

Now, let 7„ be a connected component of K n such that dist(K n ,K*) = dist(7„, By the 

maximal principle for harmonic functions, it holds that g ln > gx n for z ^ K n , and therefore, 
7* C K*. Thus, 

~ e 2 c 
dist(X n ,^) > dist( 7n ,7^) > 

by (3.7) and the definition of co- This finishes the proof of the proposition. □ 

Let A be as in Proposition C. Observe right away that A has no interior as otherwise there 
would exist A' C A with smaller logarithmic capacity which still belongs to JC. It turns out that 
<7a has a rather special structure that we describe in the following proposition which was initially 
proven in this form in [43, Theorem 3] (the method of proof in a more general form was also used 
in [20]). 

Proposition 7. Let A be as in Proposition 6. Then 



g A {z) = Re {^J 



A{r) 



where A was defined in (2.4), B is a monic polynomial of degree p — 2, and the root is chosen so 
that zy/A(z)/B(z) -t 1 as z — > 00. 

Proof. Denote by u>k the equilibrium measure of a compact set K and by I\p] the logarithmic 
energy of a compactly supported measure /x, i.e., 



I[lA = ~ 11 lo 2>\ z - T\d^{z)d^i(T). 

Then it is known that 

9a(z) = I[<*>a] + / 'log \z-t\cIua(t), 
which immediately implies that 

(3.10) 0*9A)(z) = l fd " AiT) 



2 J z-t ' 

where d z := (d x — id y )/2. Since g A = on A, it holds that 

g A (z)=Re\2^ (d z g A )(T)d T 

for any the path of integration in D. Thus, to prove (3.9), we need to show that 
/o n-v B ( z ) _ ( [ g^a(t) x 1 

for some monic polynomial B, deg(B) = p — 2. 
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Let O be a neighborhood of A. Define 

5{z): 



A{z) 



Then 6 generates a local variation of O according to the rule z i— > z := z + tS(z), where t is a 
complex parameter. Since 



(3.12) 



Iw'-z*] 



\w — z\ 



l+t 



5{w) - 5(z) 



this transformation is injective for all \t\ < to < M, where 

(3.13) M:= ma^\(6(w)-S(z))/(w-z)\<oo. 

w.zeo 

Moreover, the transformation S naturally induces variation of sets in O, E h-> E l = {z 1 : z E E}, 
and measures supported in O, \i v+ fjfi, ^(E*) = n(E). 

Let /x be a positive measure supported in O with finite logarithmic energy I\p]- Observe that 
the pull-back measure /i* satisfies the following substitution rule: d/j, t (z t ) — dfi(z). Then it follows 
from (3.12) that 



log 



l + t 



6{w) - 6{z) 



(3.14) 



-Re 



log l+t 



w — z 
5(w) - 5{z) 



dn{z)dii{w) 

d/j,(z)dfj,(w) 



for all |i| < t . Since the argument of the logarithm in (3.14) is less than 2 in modulus, it holds 
that 



(3.15) 

for all It I < to , where 



It 1 



I[fj] = -Re [tSQi) + 0{t 2 )] 
S(w) - 8{z) 



w — z 



-dfi(z)dfj,(w) 



Let now {/J-t} be a family of measures on A such that (/it)' = cja*- Then 

(3.16) <5(/i t ) -> 5(u A ) as t -> 0. 

Indeed, by the very definition of the equilibrium measure it holds that the differences I\pt] — /[wa] 
and /[w^] — J[wAt] are non-negative. Thus, 

< I[nt] - I[ua] = /[^t]-/[w A *]+/[wA*]--f[wA] 
< 7[/i { ] - J[w A t] + I[w A ] - J[w A ] 

(3.17) = Re[i(5(^)-<5(^A))+0(f2)] 

for |t| < to by (3.15). Clearly, \6(n)\ < M\fi\ by (3.13), where |/i| is the total variation of /i. Since 
Ht and ua are positive measures of unit mass, (3.17) implies that /[/it] — > I[^a] as t — > 0. The 
latter yields that /i f A u; A by the uniqueness of the equilibrium measure, 3 which immediately 
implies (3.16) by the very definition of weak* convergence. 



The measure w& is the unique probability measure that minimizes energy functional /[•] among all probability 
measures supported on A. As any weak limit point of {fit} has the same energy as oja by the Principle of Descent 
[46, Theorem 1.6.8] and (3.17), the claim follows. 
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Now, observe that a\ — for any k € {1, . . . ,p}. Hence, A* £ K, for all \t\ < to- In particular, 
this means that cp(A*) > cp(A) and therefore /[wa*] < I[^a] as cp(K) — exp{— /[wjf]}. Thus, it 
holds that 

(3.18) < J[u A ] - J[ua<] < I\pt] - I[wa<] - Re [t<5(/z t ) + (9(t 2 )] = Re [M(w A ) + o(t)] 
by (3.15) and (3.16). Clearly, (3.18) is positive only if 

(3.19) 5(u A ) = 0. 

In another connection, observe that there exists a polynomial in u, say 

B(u; z, w) — ao(z, w) + ai(z, w)u + ■ • • + a p _3(z, w)u p ~ 3 + u p ~ 2 , 
where each a^{z,w) is a polynomial in z and w, such that 

(3.20) (u> — u)A(z) — (z — u)A(w) + (z — w)A(u) = (z — w)(z — u)(w — u)B(u; z, w). 

Indeed, the left hand side of (3.20) is a polynomial of degree p in each of the variables z, w, u that 
vanishes when u — z, u = w, and z — w. Then 

M - S(w) = A(z) A(w) = A(u) 

z — w [z — u){z — w) (w — u)(z — w) 1 ' ( z ~~ u)(w — u) 

by (3.20). So, we have by the definition of 5(u>a) and (3.19) that 

B[u) := [ B { u-z, W )^)d^{ W )=A { u) ff ^(z)d. A ( w ) = / r 

J J J {z-u){w-u) \J t-uJ 

which shows the validity of (3.11) and respectively of (3.9). □ 

Having Proposition 7, we can describe the structure of a set A as it was done in [48] with the 
help of the critical trajectories of a quadratic differential [44, Section 8]. Recall that a quadratic 
differential is the expression of the form Q{z)dz 2 , where Q is a meromorphic function in some 
domain. We are interested only in the case where Q is a rational function. 

A trajectory of the quadratic differential Q{z)dz 2 is a smooth (in fact, analytic) maximal Jordan 
arc or curve such that Q(z(t))(z' (t)) 2 > for any parametrization. The zeros and poles of the 
differential are called critical points. The zeros and simple poles of the differential are called finite 
critical points (the order of the point at infinity is equal to the order of Q at infinity minus 4; for 
instance, if Q has a double zero at infinity, then Q{z)dz 2 has a double pole there). 

A trajectory is called critical if it joins two not necessarily distinct critical points and at least 
one of them is finite. A trajectory is called closed if it is a Jordan curve and is called recurrent 
if its closure has non-trivial planar Lebesgue measure (such a trajectory is not a Jordan arc or a 
curve). A differential is called closed if it has only critical and closed trajectories. 

If e is a finite critical point of order k, then there are k + 2 critical trajectories emanating from 
e under equally spaced angles. If e is a double pole and the differential has a positive residue at e, 
then there are no trajectories emanating from e and the trajectories around e are closed, that is, 
they encircle e. 

Proposition 8. Let A be as in Proposition 6 and the polynomials A, B be as in Proposition 7. 
Then (2.3) holds with [J A& being the union of the non-closed critical trajectories of the closed 
quadratic differential (—B/A)(z)dz 2 and E being the symmetric difference of the set {a\, . . . ,a p } 
and the set of those zeros of B that belong to the closure o/IJA/j. The remaining zeros of B, 
say {b\, . . . , b q }, are of even order and of total multiplicity 2(m — 1), where m is the number of 
connected components of A. Furthermore, (3.1) holds. 
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Proof. In [20, Lemma 5.2] it is shown that A is a subset of the closure of the critical trajectories 
of (— B / A)(z)dz 2 . Since (3.9) can be rewritten as 

9a(z) 

the critical trajectories of (— B / A){z)dz 2 are the level lines of gA and therefore (— B /A)(z)dz 2 is 
a closed differential. By its very nature, A has connected complement and therefore the closed 
critical trajectories do not belong to A. Since 2d z g\ = y/B/A is holomorphic in D* , all the 
non-closed critical trajectories belong to A and all the zeros of B that belong to the closed critical 
trajectories are of even order. Let us show that their total multiplicity is equal to 2(m — 1). This 
will follow from the fact that the total multiplicity of the zeros of B belonging to any connected 
component of A is equal to the the number of zeros of A belonging to the same component minus 2. 

To prove the claim, we introduce the following counting process. Given a connected compact 
set with connected complement consisting of open Jordan arcs and connecting isolated points, we 
call a connecting point outer if there is only one arc emanating from it, otherwise we call it inner. 
Assume further that there are at least 3 arcs emanating from each inner connecting point. We 
count inner connecting points according to their multiplicity which we define to be the number of 
arcs incident with the point minus 2. Suppose further that the number of outer points is p' and 
the number of inner connecting points is p' — 2 counting multiplicities. Now, form a new connected 
set in the following fashion. Fix p of the previously outer points and link each of them by Jordan 
arcs to p chosen distinct points in the complex plane in such a fashion that the new set still has 
connected complement and each of the previously outer points is connected to at least 2 newly 
chosen points. Then the new set has p' — p + p outer connecting points (p' — p of the old ones and 
p of the new ones) and p' — 2 + (p + p — 2p) = p' — p + p — 2 inner connecting points. That is, 
the difference between the outer and inner connecting points is again 2. Clearly, starting from any 
to points in E connected by a Jordan arc, one can use the previous process to recover the whole 
connected component of A containing those two points, which proves the claim. 

To prove (3.1), observe that by (3.10) and (3.11) we have that 

|^ = 2Re {n ± d z g A ) = Re 

± L \ / ±j 

where n± are the unimodular complex numbers corresponding to n±. Since the tangential deriv- 
ative of 3a is zero, so is the imaginary part of the product n±(yjB/A)±. Since, n + = — ri_ and 
{yfB]A) + = -{t/bJA)-, (3.1) follows. □ 

Propositions 6-8 are sufficient to prove Theorem 4. As an offshoot of Theorem 4 we get that the 
contour A is unique since (2.28) and (2.29) imply that all but finitely many zeros of [n/n]p converge 
to A. However, this fact can be proved directly [47, Thm. 2]. Moreover, it can be shown that 
property (3.1) uniquely characterizes A among smooth cuts making / single- valued [9, Thm. 6]. 

While Propositions 6-8 deal with the most general situation of an arbitrary finite set A, Condi- 
tion CP introduced in Section 2 is designed to rule out some degenerate cases. Namely, it possible 
for some zeros of the polynomial B to coincide with some zeros of the polynomial A. This happens, 
for example, when all the points in A are collinear. In this case, the minimal capacity cut A is 
simply the smallest line segment containing all the points in A, and the zeros of B are exactly the 
zeros of A excluding two that are the end points of A. It is also possible for the polynomial B to 
have zeros of multiplicities greater than one that belong to A. These zeros serve as endpoints to 
more than three arcs (multiplicity plus 2), see Figure 4A. However, under small perturbations of 
the set A, these zeros separate to form a set A satisfying Condition CP, see Figure 4B. 
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(A) (B) 



FIGURE 4. Both figures depict the set A for the quadratic differential 

(z 2 - b 2 )dz 2 

—. , where \a\ — 1 and b is real and depends on a. On the left-hand 

z* — {a z + a 2 -)z i + 1 

figure Arg(a) = 7r/4 which forces 6 = 0. On the right-hand figure, Arg(a) < 7r/4, in 
which case b > 0. 



4. Riemann Surface 

Let D\ be the Riemann surface of h described in Section 2.3. That is, 9t is a hyperelliptic Riemann 
surface of genus g with 2g-\-2 branch (ramification) points E and the canonical projection (covering 
map) 7r : £H — > C. We use bold letters z,w,t to denote generic points on *H and designate the 
symbol •* to stand for the conformal involution acting on the points of according to the rule 

z * = z (i-k) for Z = ^W ) fcg{o,l}. 

4.1. Abelian Differentials. For a rational function on 9\, say /, we denote by (/) the divisor of 
/, i.e., a formal symbol defined by 

(/):= £ z- £ w, 

z: /(ss)=0 w: /(w)=0O 

where each zero z (resp. pole w) appears as many times as its multiplicity. A meromorphic 
differential on *H is a differential of the form fdz, where / is a rational function on $H. The divisor 
of fdz is defined by 

(fdz) := (/) + (dz) = (/) + e- 2oo(°) - 2oo«. 

eS-E 

It is more convenient to write meromorphic differentials with the help of 

h{zM) :={-l) k (\{{z-e k )\ ' , fee {0,1}, 
^ keE ' 

where the square root is taken so h(z^) / z 9+1 — > 1 as z —> oo. Clearly, h is a rational function on 
91 with the divisor 

(n) = ^e-( 5 +l)oo(°)-( 5 +l)ooW. 

e£E 

Then arbitrary meromorphic differential can be written as fdtt, dtt :— dz/H, and respectively 

(/dO) = (/) + ( 5 -l)o (°) + ( 5 -l)oo( 1 ). 

A meromorphic differential is called holomorphic if (fdft) > (its divisor is integral). Since for 
any polynomial £ it holds that 

(£ott)= J2 ( zi0) + z(1) ) - de sW (°° (0) + 
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the holomorphic differentials arc exactly those of the form idil, deg(^) < g. Thus, there are exactly 
g linearly independent holomorphic differentials on 5H. Under the normalization (2.19), these are 
exactly the differentials <Klk- 

Let now Wi, W2 G wi 7^ W2. We denote by df2 Wl W2 the abelian differential of the third kind 
having two simple poles at Wi and W2 with respective residues 1 and —1 and normalized so 

(4.1) <f dfi Wl , W2 =0, j g {l,...,g}. 

It is also known that 



(4.2) d> dft WI ,w 2 = -2m / dSlj, j G {1, . . .,g}, 
where the path of integration lies entirely in 5H. 

4.2. Green Differential. The Green differential dG is the differential dQ.^rt modified by a 
suitable holomorphic differential to have purely imaginary periods. In fact, it holds that 

(4.3) dG(z) = h{z)dz, zem. 

Indeed, the value of the integral of dG along any cycle in £H\ |oo^ ^ U 00^} is purely imaginary 
as it is a linear combination with integer coefficients of its periods on the a- and b-cycles and 

the residues at 00^ and oo^ with purely imaginary coefficients. Thus, Re dG^j is the 

Green function for and therefore is equal to go lifted to D^°\ Hence, the claim follows from 
Proposition 7. 

For z G 9t\ {oo(°\oo (1) }, put 



(4.4) 



G(z) := f dG. 



Then G is a multi-valued analytic function on D\ \ {oq(°\ oo' 1 )} which is single-valued in 9\. 
Moreover, it easily follows from (4.3) and the fact that a\ is a branch point for *H that 

(4.5) G(z (0) ) + G(z (1) ) = (mod 2™) in D a . 

Furthermore, for any point z G Ufe=i( a fc U bfc) ^ holds that 

(- j> dG, if zGa fc \b fe , ( 2iriw k , if z G a fe \ b fc , 
T 
j> dG, if zGb fc \a fe , t 2iriT k , if zGb fc \a fc , 

where the constants u k and r k were defined in (2.17) (clearly, the integrated differential in (2.17) 
is dG). As all the periods of dG are purely imaginary, the constants oj k and r k are real. With the 
above notation, we can write 

g 

(4.7) dG = dft^D^o) + 2-Ki^TjdQ.j. 

i=i 

Indeed, the difference between dG and the right-hand side of (4.7) is a holomorphic differential with 
zero periods on a-cycles and therefore is identically zero since it should be a linear combination of 
differentials satisfying (2.19). In particular, it follows from (4.2) and (4.7) that 

,oo<°> 

(4.8) / dfi = cD + Bar. 
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Using the Green deferential dG, we can equivalently redefine $ introduced in (2.14) by 



(4.9) <&:=exp{G}, = + ... . 

Then $ is a meromorphic function on $\ with a simple pole at oo^ , a simple zero at oo^ 1 ) , otherwise 
non-vanishing and finite. Moreover, $ possesses continuous traces on both sides of each and 
that satisfy (2.16) by (4.6), and (2.15) by (4.5). Let us also mention that the pull-back function 
of $ from onto D a , which we continue to denote by <£>, is holomorphic and non- vanishing in 
D a except for a simple pole at infinity. It possesses continuous traces that satisfy 



(4.10) 



$+/$ = exp{2niuj k } on A% 
$~$+ = exp{2iTi6 k } on A fe , 



by (2.16), (2.15), and the holomorphy of $ across those cycles L k that are not the b-cycles, where 
we set 5k := Tj k if A^ = ir(hj k ), and 8k '■— otherwise. 

4.3. Jacobi Inversion Problem. Let r be a rational function on C. Then r o it is a rational 
function on with the involution-symmetric divisor, i.e., 

j j j j 

As *K is hyperelliptic, any rational function over fH with fewer or equal to g poles is necessarily of 
this form. Recall that a divisor is called principal if it is a divisor of a rational function. Thus, 
the involution-symmetric divisors are always principal. By Abel's theorem, a divisor 53j=i *j — 
Y^,j=i w j is principal if and only if fc = / and 

/ d& = (mod periods d& 



E 



In fact, it is known that given an arbitrary integral divisor Y)j—y Wj, for any vector c there exists 
an integral divisor Xw=i *j such that 



9 



(4.11) 2 / = c [mod periods 

The problem of finding a divisor X)j=i *j f° r given c is called the Jacobi inversion problem. The 
solution of this problem is unique up to a principal divisor. That is, if 

9 

(4.12) E/*- 7 — { P rmc iP a l divisor } 

is an integral divisor, then it also solves (4.11). Immediately one can see that the principle divisor 
in (4.12) should have at most g poles. As discussed before, such divisors come only from rational 
functions over C. Hence, if 2^f=i*j> a solution of (4.11), is special, that is, contains at least 
one pair of involution-symmetric points, then replacing this pair by another such pair produces a 
different solution of the same Jacobi inversion problem. However, if X^f=i *j is not special, then it 
solves (4.11) uniquely. 



STRONG ASYMPTOTICS OF NUTTALL-STAHL POLYNOMIALS 



23 



4.4. Riemann Theta Function. Theta function associated to Bq is an entire transcendental 
function of g complex variables denned by 

ri + 27riri T u>, wGC 9 . 



6 (u) := exp < ■nin T Bo,n + 2mn T u > 



As shown by Riemann, the symmetry of Bq and positive definiteness of its imaginary part ensures 
the convergence of the series for any u. It can be directly checked that 9 enjoys the following 
periodicity properties: 



(4.13) 



9 + j + Bo,fn^ = cxp / — TTirh T Bum — 2irim T u^6(itj , j, m G Z 9 



The theta function can be lifted to in the following manner. Define a vector f2 of holomorphic 
and single-valued functions in U\ by 

(4.14) fi(z) := j dn, z G VK. 

J a± 

This vector-function has continuous traces on each side of the a- and b-cycles that satisfy 

,. lr , x+ f -Bnek on a fe , 

(4.15) ={ _ k€{l,...,g}, 



Sk on b 



A- 



by (2.19) and (2.20). It readily follows from (4.15) that each fi^ is, in fact, holomorphic in 9t\ hj~. 
It is known that 



s- 1 , 

(4.16) 0(u)=O <^ u = 22^1 (tj) + K (mod periods dfi 



for some divisor X)f=i where K is the vector of Riemann constants defined by (K)j := ((Bn) 
l)/2-Zk*ifa*to7<Mk,je{l,... t g}. 

Let y]j—i tj and Y)j—i Zj be non-special divisors. Set 

4.17 9 z;Vt 3 ,yz 3 :=-^ = ~ 
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It follows from (4.15) that this is a meromorphic and single-valued function in SH (multiplica- 
tively multi- valued in Furthermore, by (4.16) it has a pole at each Zj and a zero at each tj 
(coincidental points mean increased multiplicity), and by (4.13) it satisfies 

(4.18) 9+ (z;^t„^z,) = 0- (z;^t^z/) cxp j 2 ™E ~ <M%)) 

for zeafeUUZj-IJUtj}. 

If the divisor Ylj=i^j ( r esp. J2j=i z j) m (4-17) is special, then the numerator (resp. denomi- 
nator) is identically zero by (4.16). This difficulty can be circumvented in the following way. Let 
Wi , W2 G SH \ {w} for some w G 5H. Set 

9 (Yz(z) - fi(wi) -{g- l)n (w*) - K 

(4.19) 6(z;w 1 ,w 2 ):- V 



(n(z) ~ n(w 2 ) -(g- i)n (w*) - k 
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Since the divisors w_,- + (<? — l)w* are non-special, (■;wi,W2) is a multiplicatively multi-valued 
meromorphic function on U\ with a simple zero at Wj, a simple pole at w 2 , and otherwise non- 
vanishing and finite. Moreover, it is meromorphic and single-valued in d\ and 

(4.20) 6+ (z; wi.wj) = 0- (z; w x , w 2 ) exp (2^(fi fc (w 2 ) - fi^Wi)) } 

for z G aj; \ {wi,W2}. Observe that the jump does not depend on f2(w*). Hence, analytic 
continuation argument and (4.20) immediately show that 0(-;wi,W2) can be defined (up to a 
multiplicative constant) using any divisor Ef=i as l° n S as w * + E?=i *j i s non-special and that 



©O^E^-E^) i r 6(z;tj 



n 



© (w; E tj , E z j ) ^ Q (w; t j , zj ) 

for any w fixed and satisfying {w} n {UZj Jutj} = 0. Let us point out that even though the 
construction (4.17) is simpler, it requires only non-special divisors, while this restriction is not 
needed for (4.19). 

5. Boundary Value Problems on L 

This is a technical section needed to prove Proposition 3. The results of this sections will be 
applied to logarithm of p £ Wa, which is holomorphic across each arc comprising A. However, 
here we treat more general Holder continuous densities as this generalization comes at no cost 
(analyticity of the weight p will be essential for the Riemann-Hilbert analysis carried in Sections 7- 
8). In what follows, we describe properties of 

(5.1) := I ^tffl.,,., z£K\L, 

4ttz J l 

for a given function ip on L. Before we proceed, let us derive an explicit expression for d£t z z - . To 
this end. set 



(5.2) 



H k (z):=d> =2/ k€{l,...,g}. 

Ja t W" Z J A a W- Z h[w) 



Clearly, each H k is a holomorphic function on £H \ that satisfies 
(5.3) HH+ - HHfr = Am on a fc 

due to Sokhotski-Plemelj formulae [24] as apparent from the second integral representation in (5.2). 
Then, using functions H k , we can write 



(5.4) dfW(w) = -^-dn(w) - Y(hH k )(z)<m k (w). 

w — z * — ' 

k=l 

5.1. Holder Continuous Densities. Let -0 be a function on L \ E with Holder continuous 
extension to each cycle L k and ^ be given by (5.1). The differential dVt z z * plays a role of the 
Cauchy kernel on with a discontinuity. Indeed, it follows from (5.4) that 

j j — i j — i 

Each function is holomorphic in £R\ (Lj U {oo 1 - ', oo' 1 )}) with Holder continuous traces on 
Lj that satisfy 

(5.6) *+ - *7 - 0. 
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Clearly, (e) = for e g E\L r Moreover, it holds by (5.6) and the identity {z)+$ Lj (z*) = 
that 

(5.7) ^( z W)^i^L^.( e ) as 

for univalent ends e E E P\ Lj. To describe the behavior of Vf^ . near trivalent ends, recall that A 
splits any disk centered at e of small enough radius into three sectors. Two of these sectors contain 
part of Aj in their boundary and one sector does not. Recall further that H(z) changes sign after 
crossing each of the subarcs of A. Thus, it holds for trivalent ends e g E n Lj that 

(5.8) y Lj ( z W)^±t}Lty Lj ( e ) as z e, 

where + sign corresponds to the approach within the sectors partially bounded by Aj and the 
— sign corresponds to the approach within the sector which does not contain Aj as part of its 
boundary. 

Similarly, each ^J ak is a holomorphic function in 9t\ (&k U {oo*- -*, oo^}) with Holder continuous 
traces on a& that satisfy 



(5.9) ¥+-¥-=0Vdn* 



by (5.3). Analogously to (5.7), one can verify that \I/ afc (e) = for e g E \ a^ and 
(5.10) *a fc (^ (fc) ) -> I ipdVk as z^eesi k r\E. 



(5.12) = 



* JL 

Combining all the above, we get that ^ is a holomorphic function in 91 \ L including at oo^ ^ 
and oo W where it holds that 

(5.11) *(oo (0) ) = -#(oo (1) ) = -t t I ipdti - — I ipdG 

2 J L 4tti J l 

by (4.7). Moreover, it has Holder continuous traces on both sides of (L U (J afe) \ E that satisfy 

tp, on L\E, 

-§ L ipdVt k , on a fe \_E, 

according to (5.6) and (5.9). Finally, the behavior at e g E can be deduced from (5.7), (5.8), and 
(5.10). 

5.2. Logarithmic Discontinuities. Assume now that ip has logarithmic singularities at e g E, 
which, obviously, violates the condition of global Holder continuity of tp on the cycles L k . However, 
global Holder continuity is not necessary for \I> to be well-defined. In fact, it is known that the 
traces 1 $> ± are Holder continuous at t g L as long as tp is locally Holder continuous around this 
point. Thus, we only need to describe the behavior of VP near those e g E where tp has a singularity. 

Let a be a fixed univalent end of A and Aj be the arc incident with a. Further, let ip be a fixed 
determination of alog(- — a) holomorphic around each arc A^ (except at a when k — j), where 
a is a constant. As before, define \& by (5.1). It clearly follows from (5.5) that we only need to 
describe the behavior of around a as the behavior of the other terms is unchanged. To this 
end, it can be readily verified that 

(5.13) <M-)" nW f m dt 



2ni J A . t-zh+{t)' 
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Denote by U a ,s a ball centered at a of radius d chosen small enough that the intersection AjDU^s is 
an analytic arc. Denote also by U^g the maximal open subset of U a ^\Aj in which ip is holomorphic 
and Aj C dU^g if Aj is orienter towards a and A^ C dU^ s if Aj is oriented away from a. Set 

(5.14) arg(a — z) := arg(z — a) ± 7r in U^ s . 

It can be readily verified that thus defined log(a — z) := log \a — z\ + i arg(a — z) is holomorphic 
in U a .s \ Aj. Then, arguing as in [24, Equations (8.34)-(8.35)], that is, by identifying a function 
with the same jump across Aj as the one of integral in (5.13), we get that 

1 f 4>(t) dt a \og(a — z) r , , . , . . 

(5.15) / , . , = — -7 1- \ terms that are holomorphic at a \ 

K ' 2m J A t- zh+(t) 2 h(z) 1 ' 

in U a> s \ Aj. Multiplying both sides of (5.15) by H, we get that (5.7) is replaced by 

(5.16) ^ L] {z (k) ) ~ (-l) fe |log(a- z) -> as z -> a. 

Let now b be a trivalent end of A and ip be a fixed determination of alog(- — b) analytic across 
A, where a as before is a constant. Further, let A^j be the arcs incident with b. Fix I € {1, 2, 3} 
and let log(6 — z) be defined by (5.14) with respect to Atj. Then (5.15) still takes place within 
the sectors delimited by A^/ U Abj+i and Abj U A& where I ± 1 is understood cyclicly within 
{1, 2, 3}. However, within the sector delimited by A^j+i U A& the right-hand side of (5.15) has 
to be multiplied by —1 to ensure analyticity across A^z+i U Ab t i-±. Then, multiplying both sides 
of (5.15) by ft, we get that 

(5.17) *L t ,(z (fc) )T(-l) fc | log(6-z)->0 as z -> b, 

where — sign corresponds to the approach within the sectors delimited by A^ U A^j+i and A;,.; U 
Ab : i-i, and the — sign corresponds to the approach within the sector delimited by the pair A^+i U 
Ab i—i- Hence, the behavior of ^ near b is completely determined by the behavior of the sum 

£?=!**».,• 

5.3. Auxiliary Functions. For an arbitrary f £l 9 set xpg to be a function on L such that 

2m(x), on bfc, 



Define further 



otherwise. 



(•>.Im .S', ( z) : =' >X P|^ 4 r,-r/<) zz \. z^!',L 



Then is a holomorphic and non- vanishing function in IK, with Holder continuous non- vanishing 
traces on both sides of each a- and b-cycle that satisfy 



exp {27r«(a?) fe } on b 



k: 



(5-19) St = SZ 

y exp { - 2m {Bnx) k ) on a fc , 

for k £ {1, . .. ,g} by (5.12) and (2.20). Observe also that 

(5.20) %=n4f= 



fe=l 



where, as before, {efe}| =1 is the standard basis in 
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Now, let p £ Wa and log p be a fixed branch holomorphic across each A& in A. Define 

1 



(5.21) -Siogp(z) := exp i — <b log(po Tr)dn z ^ \ , zeSH\L 



L 



Then, as in the case of (5.18), it follows from (5.12) that Siogp is a holomorphic function in 9\\ L 
and 



(5-22) S+ gp = S^ p 



po7r, on L\E, 

ex P{ ~ § L l °g(P 07r ) dn k}, on a k \E, ke{l,...,g}. 



If a is a univalent end of A and Aj is the arc incident with a, then p\A j { z ) — w j( z )( z ~ a ) Qa : 
where Wj is holomorphic and non- vanishing in some neighborhood of Aj and (z — a)" a is a branch 
holomorphic around Aj \ {a}. Let (a — z) Qa be the branch defined by (5.14). Then (a — z) tta = 
(z — a) Qa exp{±a a 7ri} in U^ s , where the latter were defined right before (5.14). Thus, it holds by 
(5.7) and (5.16) that 



(5.23) S^ ORP {z)/p{z) -» exp |±a g 7ri - e &k {a, z) j> log(p o 7r)dOfc j as U^ s 3 z ^ a, 

where e a( .(a,z) = if a ^ a^, and e afc (a,z) = ±1 if a £ aj. and z — > a G a^ . If & is a trivalent 
end of A, let Abj, j £ {1, 2, 3}, be the arcs incident with &. Further, let Sbj be a sector delimited 
Abj±± within a disk centered at b of small enough radius, where j '± 1 is understood cyclicly within 
{1, 2,3}. Then it follows from (5.8) that 

(5-24) ^ P (^)^ ' |A6 ' J - lW T,r i(b) as 5 b , 9 ^6. 

Finally, we can deduce the behavior of S\ osp at oo^ -* from (5.11) in a straightforward fashion. 

Now, let logh + be a fixed branch continuous on A \ E. Define S\ og h as in (5.21) with logp 
replaced by \ogh + . Then S\ og h enjoys the same properties S\ ogp does except for (5.23) and (5.24) 
as h is not holomorphic across A \ E unlike p. In particular, it holds that 

h + o 7r, on L\E, 

exp { - § L log(h + o ir)dQ k } , on a k \E 7 k G {1, . . . ,g}. 

Moreover, it can be easily verified that (5.23) gets replaced by 



(5.25) S+=S lo , 



(5.26) Sf ogh (z)/h(z) -> exp j^Ty - ^£* k { a i z )j> log(h, + o 7r)dfi fc ^ as 



with — sign used when Aj is oriented towards a (one needs to take UT'g = C/ a .,5 \ Aj and U s = 
in (5.16)) and + sign used when Aj is oriented away from a (U~ s = U a ,s \ Aj and U^g = )- 
Similarly, (5.24) is replaced by 

(5.27) S£ gfc //i->exp{±7ri/2} as z -> 6, 

where the sign + corresponds to the case when the arcs incident with & are oriented towards b and 
the sign — corresponds to the other case (this conclusion is deduced from (5.8) and (5.17) applied 
to each arc incident with b). 



6. Szego Functions 
In this section we prove Propositions 1, 2, and 3. 
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6.1. Proof of Proposition 1. It follows from the discussion in Section 4.3 that a solution of 
(2.22) is either unique or special; and in the latter case any pair of involution-symmetric points 
can be replaced by any other such pair. According to the convention adopted in Definition 2, 
we denote by X^f=i ^n,j the divisor that either uniquely solves (2.22) or solves (2.22) and all the 
involution-symmetric pairs are taken to be oo^ -* -I- oo^ 1 -*. 

Let Ylj=i "j be as in (2.13) and lu,t,c p be as in (2.21). Notice that by (4.8) it holds that 

9 / >oo(0) / — *\ 

(6.1) J dti + i J d& = c p + (n + i) (p + Bqt) (mod periods d(l J . 

Then if 

g g — l 

J2 t W j = + fcoo (0) + (I - k)oo W , 

i=i i=i 

where {tj}^Z[ C 5K\ {ooW.ooW}, it holds by (6.1) that 

a g-i 

J2 = J2 *j + ( fc + *)°° (0) + (* ~ fc ~ 

i=i i=i 

for each z € {— fc, — fc}. The uniqueness of the solutions for i — —k, I — k immediately follows 
from the fact that 53f=i ^n+i,j is n °t special for these indices. 

Now, let Y^j=i^"n,j be the unique solution of (2.22) that does not contain oo^ fe \ k 6 {0, 1}. If 
Sf=i ^n— (— w ere not the unique solution, it would contain at least one pair c©W + oo^ and 
therefore Ylj=i ^n,j would contain oo^ fc ^ by the first part of the proof. Thus, X)f=i t ri -(-i) fc .j solves 
(2.22) uniquely, and it only remains to show that 

(6-2) {W}- =1 n{t„-(-i )tj }; =1 = 0. 

Assume the contrary. For dcfinitcncss, let g' < g be the number of distinct points in the divisors 
Y]j—i t n ,j and Sf=i t n — (_ and label the common points by indices ranging from g' + 1 to g. 
If (6.2) were false, then it would follow from (2.22) and (4.8) that 



a rt n ,.oo<°> 

E 

3=1 



[ dtt- (-l) fe / dti = (mod periods dti) . 



That is, the divisor 



£t nJ - ^t„_ ( _ 1)fcj - (-1)W°> + (-l)*oc« 

would be principal. However, since g' + 1 < <?, such divisors come solely from rational functions 
over C and their zeros as well as poles appear in involution-symmetric pairs. Hence, the divisor 
Ef=i would contain an involution- symmetric pair or oo^ k K As both conclusions are impossible, 
(6.2) indeed takes place. This completes the proof of Proposition 1. 

6.2. Proof of Proposition 2. Let N" C N' be such a subsequence that the divisors J2i=i ^n+j,i 
converge to a divisor $2?=i w « as ^" 3 « — > oo for a fixed index j € {— Iq — k, . . . , h + Then 
the continuity of implies that 
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where t?& := Yli=i ) (recall also the convention that all the paths of integration belong to 9t 

and therefore the right-hand side of the equality above does not depend on the labeling of ^ ^n+j.i 
and Y^b\ )■ Hence, it holds that 

9 

lim (c p + (n+j)(uj + But)) = (w,) - v b , 

i=l 

where, from now on, all the equivalences are understood mod periods d£l. Set / := lo — 1\. Assume 
first that I > 0. Then, analogously to the previous computation, we have that 

g— 2k— lo— h 

lim {c p + n{uj + BuT))= E ^ + ^ (°° (0) ) - v b 
since ft («< )) = -ft 

In what follows, we assume that l + 2j > 0, otherwise, if l + 2j < 0, each occurrence of oo^ and 
I + 2j needs to be replaced by oo^ and —(1 + 2j), respectively. Then (4.8) and the just mentioned 
anti-symmetry of ft yield that 

g-2fc-Z -(i 

lim (c p + (n + j)(cj + B n f))= V CL(t z ) + (l + 2j)n(^°A -v b . 

1=1 

Hence, it is true that 

g g—2k — lQ — li 

^!J(w,)e ft(t;) + (Z + 2j)ft(oo (0 )). 

i=l i=l 

Therefore, for any collection {u,}i=i' £ ^ C C it holds by Abel's theorem that the divisor 

g—2k — lo—li li+k—j g 

E *+ E (-f ) +^ 1) ) + ( /+2 ^ (0) -E^ 

2=1 2 = 1 2 = 1 

is principal (?i + fc — j needs to be replaced by + fc + j when I + 2j < 0). As the integral 
part of this divisor has at most g elements, the divisor should be involution-symmetric. However, 
if ^^~ 1 2fc_ ' 0_il tj + (I + 2j)oo(°) is non-void, it is non-special, and therefore J2i=i w * i s equal 

to J2W~ lo ~ h U + Eti* _i («! 0) + + + 2j)oo (0) ; if it is void, £f =1 w, is an arbitrary 
involution-symmetric divisor. In any case, this is exactly what is claimed by the proposition. 
Clearly, the case I < can be treated similarly. 

To prove the last assertion of the proposition, observe that the divisors X)f=i ano - Tlj—i w i 
are connected by the relation 

a _ <.oo (0) a 

E^(t,)-(-i) fc / dfi = x;6(wi). 

Hence, by Abel's theorem the divisor Y^j=i^j ~ Ef=i w i — ( — l) fe oo' ^ + (— l) fe oo' 1 ^ is principal. 
Since Ylj=i *j + oo( 1 ~ fe * 1 is non-special, the claim follows as in the end of the proof of Proposition 1. 

6.3. Proof of Proposition 3. Any vector u € C 9 can be uniquely and continuously written as 
x + Bny, x , y € M. 9 , since the imaginary part of Bn is positive definite. Hence, we can define 

(6.3) x n +B n y n ~ E L' ^ = E (^(W) - " W 5 )) ■ 

2=1 ^ fc 5 1) 2=1 
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As the image of the closure of *H under is bounded in C 9 , it holds that 
(6.4) \x n \, \y n \ < const, 

independently of n, where |c| 2 := Y^k=i \{^)k\ 2 - Set further 



c 



p- 



Then it follows from the very choice of J2j=i^n.j, see (2.22), that there exist unique vectors 

jn > ITT-n G % 9 Such that 

(6.5) x p + nti = $ n +j n and y p + nf = y n + m n . 
Therefore, we immediately deduce from (6.4) that 

(6.6) |m n — m„_i|, |(2n — l)r — m n — m„_i|, | nr — m n | < const, 
independently of n. 

Let now S? and be defined by (5.18). Then it is an easy consequence of (6.6) and (5.20) 
that 

(6.7) < const. < ISrAjSrfi^ I , ISl^/S^S^^ | , \S 7% JS?\ < const. < oo 
uniformly in 5H. Notice also that 

f exp{-27rmT fe } on b fc , 

(6.8) (S, Tl JS?) + = (S rn JS?r \ 

{ cxp { - 2m (B n (m n - nf)) k } on a fc . 

Using definitions (4.17) and (4.19), set 

g 

e n (z) : = e (^XXj-X^r) and 6n(z) : = LI 9 ( z;t »j' 5 . 



j 

3 = 1 



where the first formula is used for non-special divisors ^n,j and the second one otherwise. Then 
n is a meromorphic function in *H with poles at of , zeros at t nj (as usual, coincidental points 
mean increased multiplicity), and otherwise non-vanishing and finite. It also follows from (4.20), 
(4.18), and (6.3) that 

(6.9) 9+ = 9~ exp{-2?ri (x n + B n y n ) k } on a fc . 

Finally, let Si og h and 5iogp be defined as in Section 5.3 with the branch of the difference 
logp — \ogh + chosen to match the one used in (2.21) to define c p . Set 

(6.10) Sn :=^^Q n . 

Then S n is a meromorphic function in *H\ L and S n <& n is meromorphic in JK\ L by (2.16), (5.22), 
(5.25), (6.8), (6.9), and (6.5). Clearly, the same equations also yield that S n $ n satisfies (2.25). 
Finally, (2.26) follows from (5.23) and (5.24), (5.26) and (5.27), reciprocal symmetry of S\ ogp and 
S\ og h on different sheets of Cft, and the properties of 0„ . 

Now, let S be as described in the statement of Proposition 3. Then by the principle of analytic 
continuation S/S n is a rational function over 1H with the divisor — 53f=i - Since rational 
functions have as many zeros as poles, the divisor ^ tj has exactly g elements. Further, as 
explained in Section 4.3, the principal divisors with strictly fewer than g + 1 poles are necessarily 
involution-symmetric; that is, they come from the lifts of rational functions on C to *H. It also 
follows from Proposition 1 that Y)j—x t n j consists of a non-special part and a number of pairs 
ocA 1 ) + oo(°). Hence, Yl^j nas the same non-special part as Y^=i^n,j an d the same number of 
involution-symmetric pairs of elements. Due to Proposition 1, the latter means that J^tj solves 



STRONG ASYMPTOTICS OF NUTTALL-STAHL POLYNOMIALS 



31 



(2.22). Lastly, as all the poles of the rational function S/S n are equally split between oo*- ' and 
oo^ 1 ', this is a polynomial. 

It remains to show the validity of (2.27). It follows from the definition of S n and (6.7) that we 
only need to estimate 

e n _!(z) e n (ooW) 
e„(z) e^ooW)' 

To this end, denote by £° and €.\ the closures of { J^t^j} N and {Etn-i,j} neN m * nc 
D\ 9 /Y, g -topology. Neither of these sets contains special divisors. Indeed, both sequences consists 
of non-special divisors and therefore we need to consider only the limiting ones. The limit points 
belonging to £° are necessarily of the form 

g-2k-l k 

E t i+ E(4 0) +4 1) )+^ (1) > 

i=l i=l 

where E?=i ~ KO^OI < °°j i s non-special and {zi}^ =1 C C. If k > 0, Proposition 2, applied 
with Iq = 0, h = 1, and j = — 1, would imply that contains divisors of the form 

g-2k—l k' 

E U + (4 0) + +(k-k'- l)oo<°> + (I + 1 + k - fcV (1) 

i=l i=l 

< k' < k — 1. In particular, it would be true that I + 1 + k — k' > 2, which is impossible by the 
very definition of N £ . Since the set $} e can be examined similarly, the claim follows. 
Hence, given 

G £e: we can define ©(z; E tj>E via ( 417 )- % the 

very definition of 

£g, it holds that 



< 



e(°o w ;5>.£*S i: 



< 00. 



Moreover, compactness of £ £ and the continuity of E "(tj) with respect to E *j imply that there 
are uniform constants c(£*) and C(£g) such that 



for any E*j Analogously, observe that the absolute value of 8„_i/0„ is bounded above in 

U\ n ,e as it is a meromorphic function in with poles given by the divisor ^t n ,,, The fact that 
this bound is uniform follows again from continuity of f2 and compactness of . 

For future reference, let us point out that a slight modification of the above considerations and 
(6.7) lead to the estimates 

(6.11) \SJS n (cX>)\ , |/fcS£_i/S^i(oo)| < C e , e ,p < OO 

that holds uniformly in D* \ U e£ £;{z : \z — e\ < e} for all n £ N e . 

7. RlEMANN-HlLBERT PROBLEM 

. f ct> m \ 

In what follows, we adopt the notation <j) m<J3 for the diagonal matrix _ , _ TO , where 03 



<j>- 

is the Pauli matrix <t 3 = ( ^ ^1 ) ' M° reover i f° r brevity, we put 7a := cp(A). 
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7.1. Initial Riemann-Hilbert Problem. Let ^ be a 2 x 2 matrix function. Consider the fol- 
lowing Riemann-Hilbert problem for y (RHP-J^): 

(a) y is analytic in C \ A and lim y(z)z~ n<J3 = X, where I is the identity matrix; 



(b) y has continuous traces on each A& that satisfy y+ = J 7 - 



1 P 
1 



(c) y is bounded near each e £ E \ A and the behavior of y near each e € A is described by 
( ° ! 1 jz-ef 1 ' 1[ " ' 

O ! } J 08 ! 2-6 ! I . if ... = (i. as IV i : :> ,-. 



1 be 



z-e 



o ( \ \ J , if « e >0, 

The connection between RHP-} 7 and polynomials orthogonal with respect to p was first realized 
by Fokas, Its, and Kitaev [22, 23] and lies in the following. 

Lemma 1. If a solution of RHP-} 7 exists then it is unique. Moreover, in this case deg(g„) = n, 
R n -i{z) ~ z~ n as z oo, and the solution of RHP-} 1 is given by 



(7.i) y 



Qn Rn 

mn-lRn-l 



where m n is a constant such that m„_ 1 i? n _ 1 (z) = z "[1 + o(l)] near infinity. Conversely, if 
deg(<7„) = n and R n -\{z) ~ z~ n as z — > oo, then y defined in (7.1) solves RHP-} 7 . 

Proof. In the case when A = [—1,1] and p > on A this lemma has been proven in [32, Lemma 2.3]. 
It has been explained in [12] that the lemma translates without change to the case of a general 
closed analytic arc and a general analytic non- vanishing weight p, and yields the uniqueness of the 
solution of RHP-} 7 whenever the latter exists. For a general contour A the claim follows from the 
fact that R n = R n k, where 

D / x f q n {t)p{t) dt f q n {t)w k (t)(t-a) a *(t-b) a » dt 



t — z 2tti /a t — z 2ni 

and therefore the behavior of R n near e € A is deduced from the behavior there. On the other 
hand, for each arc A ej - incident with e £ E\A (see notation in (2.8)), the respective function R n k 
behaves as [24, Section 8.1] 



Pe,j(e) 

2,7, " ' > 



log(z - e) + R* eJ (z), 



where the function R* e ■ has a definite limit at e and the logarithm is holomorphic outside of A e j . 
Using (2.10), we get that 

R(z) = p^e) ^ _ e) + p^e) ^ _ e) + p^e) ^ ^ _ &) + ^ 

where R* has a definite limit at e and arg e ■ (z — e) has the branch cut along A e j . Thus, f} 7 is 
bounded in the vicinity of each e € E \ A. 

Suppose now that the solution, say y = [}jfc]| fc=1 , of RHP-} 1 exists. Then 3^ii = z n + lower 
order terms by the normalization in RHP-J^a). Moreover, by RHP-y(b), 3^n has no jump on A 
and hence is holomorphic in the whole complex plane. Thus, J^ii is necessarily a polynomial of 
degree n by Liouville's theorem. Further, since 3^12 = 0(z~ n ~ 1 ) and satisfies RHP-J 7 (b), it holds 
that 3^12 is the Cauchy transform of y±ip. From the latter, we easily deduce that ^ii satisfies 
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orthogonality relations (1.10). Applying the same arguments to the second row of y, we obtain 
that 3^21 == Qn-i and 3^22 = w n _i-R n _i with m„_i well-defined. 

Conversely, let deg(q„) = n and R n -i(z) = 0(z~ n ) as z — >• oo. Then it can be easily checked 
by the direct examination of RHP-^(a)-(c) that y, given by (7.1), solves RHP-J^. □ 

7.2. Renormalized Riemann-Hilbert Problem. Suppose now that RHP-^V is solvable and y 
is the solution. Define 

(7.2) T :=yX n<T3 y®~ na3 , 

where, as before, we use the same symbol $ for the pull-back function of $ from D^°\ By (4.9) it 
holds that fim z _ s . 00 z/$(z) — 7a and therefore 

(7.3) lim T(z) = lim lA n,T3 y(z)z- n ' 73 (z/$(z)) n(T3 =1. 
Moreover, it holds by (4.10) that 

(7.4) (<jj+)~ no "3 = ^([)-^-" '3 e -27rinw fc (T3 

on each A^. Finally, on each A/, we have that 



1 o i ; v o ($+/$ 



— \n 



(7.5) 



e 2irinSk ($ _ ) _2 r. 

where the second equality holds again by (4.10) and are defined right after (4.10). Combining 
(7.3) — (7.5), we see that T solves the following Riemann-Hilbert problem (RHP-T): 

(a) T is analytic in D a and T(oo) = X; 

(b) T has continuous traces on (J A& U (J A^ that satisfy 

-2mnu k a 3 Qn each A a_ 

7+ = T— \ / p 2irin8 k ( fo + \-2n „2mnSk 



2-KinSh (fi)+\— 2n e 2Trin8 k p 

^^(j-j-an ) oneach A *5 



(c) T has the behavior near each e £ £ as described in RHP-3^(c) only with D* replaced by 
D a . 

Trivially, the following lemma holds. 

Lemma 2. RHP-T is solvable if and only if RHP-J^ is solvable. When solutions of RHP-T and 
RHP-} 7 exist, they are unique and connected by (7.2). 

7.3. Opening of Lenses. As is standard in the Riemann-Hilbert approach, the second transfor- 
mation of RHP-} 7 is based on the following factorization of the jump matrix (7.5) in RHP-T(b): 
1 \ / e ^in5 kp \ / x o 

($-)" 2 »/p 1 J y - e -^in8 k j p o j y ($+)"2n/ p 1 

where we used (4.10). This factorization allows us to consider a Riemann-Hilbert problem with 
jumps on a lens-shaped contour E (see the right-hand part of Figure 5), which is defined as follows. 
For each trivalent end e, let A. e j, j £ {1, 2, 3}, be the arcs in A incident with e. For definiteness, 
assume that they are ordered counter-clockwise; that is, when encircling e in the counter-clockwise 
direction we first encounter A ej i, then A ej 2, and then A e> 3. Assume also that e > is small 
enough so that the intersection of the disk centered at e of radius e, say U e<e , with any A ej is a 
Jordan arc and the disk itself is contained in the domain of holomorphy of each w e j (see (2.11)). 
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FIGURE 5. The left figure: the arcs A ej introduced in the construction of the lens 
E near a trivalent end. The right figure: the full lens E consisting of the arcs Afc(not 
labeled), the arcs A%, and the outer arcs A^± (the choice of ± is determined by the 
chosen orientation of the corresponding arc Afe), and the domains Qk± (shaded areas, 
not labeled). 



Firstly, let A ej , j S {1,2,3}, be three open analytic arcs incident with e and some points on the 
circumference of t/ e , e placed so that the arc A ej splits the sector formed by A e:) _i and A ej+ i, 
where we understand j ± 1 cyclicly within the set {1, 2, 3}. We orient the arcs A ej so that all the 
arcs including A ej are simultaneously oriented either towards e or away from e (see the left-hand 
side of Figure 5). Secondly, let A& be an arc with one univalent and one trivalent endpoint, say e. 
Then we chose open analytic arcs Ak± C D a so that Aj;UA|i + UA e j and AfcUAfc_UA e) z (I = j + 1 
if Afc is oriented towards e and / = j — 1 otherwise) delimit two simply connected domains, say 
and fife-, that lie to the left and right of Afe (see the right-hand side of Figure 5). We oriented 
Afe± the way Afe is oriented and assume that they lie within the domain of holomorphy of w^. The 
cases where Afe is incident with two univalent ends or two trivalent ends, we treat similarly with 
the obvious modifications. Finally, we require all the arcs Afe, Afe-t, A^, and A ej - to be mutually 
disjoint, in particular, we have that A^ n fij± = for all possible pairs of k and j. 
Suppose now that RHP-T is solvable and T is the solution. Define S on C \ E by 



(7.6) 




Then ,5? solves the following Ricmann-Hilbert problem (RHP-5): 

(a) S is analytic in C \ £ and S(oo) = X; 

(b) S has continuous traces on £ that satisfy 

(1) S+ = S_ e - 27rmuJka3 on each A£; 
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1 \z-e\ 
log|z-e| log|z-e| 
logjz-ej logjz-ej 
1 1 



if 


a e 


<o, 


as C 




if 


a e 


= 0, 


as { 




if 


a e 


>o, 


as z - 


-> e outside the lens E, 


if 


Oie 


> o, 


as z - 


e inside the lens S. 



(7.7) 3T*{ Z ) 



(c) S is bounded near each e 6 E \ A and the behavior of S near each e € A is described by 
, 1 \z-e\ 

O 

°> 1 1 

/ |«-e|-°" 1 
U ^ |z-e|- Q -= 1 

Then the following lemma holds. 

Lemma 3. RHP -S is solvable if and only if RHP-T is solvable. When solutions of RHP -S and 
RHP-T exist, they are unique and connected by (7.6). 

Proof. By construction, the solution of RHP-T yields a solution of RHP-5. Conversely, let 5f* 
be a solution of RHP-5. It can be readily verified that 3F* , obtained from Sf* by inverting (7.6), 
satisfies RHP-T(a)-(b). Denote by 3~* k the jfc-entry of ST* , j, k g {1, 2}. The appropriate behavior 
of £?* 2 near the points of E follows immediately from RHP-6>(c) and (7.6). Thus, we only need to 
show that Sff\ — 0(1) in the vicinity of E and only for e € A. Observe that by simply inverting 
transformation (7.6), we get that 

if a e < 
if a e — 0, 

if a e > and z is inside the lens, 
if a e > and z is outside the lens, 
for j = 1, 2. However, each 2?* x solves the following scalar boundary value problem: 

f e -2-*** on UA^, 

1 j 9 9 \ e 2 ™«*($ + )- J " on UAfe, 

where is a holomorphic function in D a . It can be easily checked using (4.10) that $~~™ is 
the canonical solution of (7.8). Hence, the functions <j>j := ^*$™, j = 1,2, are analytic in 
C \ E. Moreover, according to (7.7), the singularities of these functions at the points e € E 
cannot be essential, thus, they are either removable or polar. In fact, since <j>j(z) = 0(1) or 
<f>j(z) = 0(\og\z — e|) when z approaches e outside of the lens E, cf>j can have only removable 
singularities at these points. Hence, (f)j(z) — 0(1) and subsequently ^* = 0(1) near each e e E 
(clearly, these functions have the form where q is any polynomials of degree at most n). □ 

8. Asymptotic Analysis 

8.1. Analysis in the Bulk. As <3? -2n converges to zero geometrically fast away from A, the 
second jump matrix in RHP-5(b) is close to the identity on |J A^-i-. Thus, the main term of the 
asymptotics for in D a is determined by the following Riemann-Hilbert problem (RHP-AA): 

(a) Af is analytic in D a and Af(oo) = X; 

(b) Af has continuous traces on (J A k U \J that satisfy RHP-5(bl)-(b2). 

As usual, we denote by S n and «S* the pull-back functions of S n on £H defined in Proposition 3. 

Lemma 4. IfnG N e , then RHP-A/" is solvable and the solution is given by 

l/S n (oo) \ / S n hS; t 
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Moreover, det(yV) = 1 on C, and it holds that Af behaves like 
D a 9 z — > e /or univalent and trivalent ends e £ E, respectively. 



|z-e|- 1/4 



2 



Proof. Observe that whenever n £ N e it holds that S n (oo)S^ l _ 1 (oo) ^ by the construction and 
therefore Af is well-defined for such indices. Since S n and h^S'^ are holomorphic function in D a 
by (2.26), Af is an analytic matrix function there. The normalization 7V(oo) = X follows from the 
analyticity of S n and £>* at infinity and the fact that (/i$)(z) = 1/ja + 0(1/ z). Further, for any 



A£ we have that 



l/5„(oo) \ / £- e - 2 ™ /z(S* 

7a/3U(oo) / V (5„_ 1 /$)- e - 2 ™^ /i($5; 



where we used (2.16), analyticity of S n $ n across the a-cycles, and the fact that 5+ / S~ 
(S*)-/(S*)+. Moreover, for each A k it holds that 

5± = (S* n rexp{~2irinS k }(h + /p) 

by (4.10), (2.25), and (2.15). Then 

l/5 B (oo) \ / -(hS* n )-e- 2mn5k /p S-e 27rinS «p 

7a/3U(oo) ) [ -(hS*_^)-e- 2 ^/p (5 n _ 1 /$)"e 2 ""' 5 V 

_ e -2-Kin5 h l p o 



Af + = 

= Af- 



on Afc, again by (4.10), (2.25), (2.15), and since h~ = —h + there. Thus, Af as defined in (8.1) 
does solve RHP-TV. Equations (8.2) readily follow from (2.26). Finally, as the determinants of the 
jump matrices in RHP-TV(b) are equal to 1, det(Af) is a holomorphic function in C \ E. However, 
it follows from (8.2) that 

det(Af)(z) < const. \z - e|~ 1/2 as z-)ee£. 

Thus, det(jV) is a function holomorphic in the entire extended complex plane and therefore is a 
constant. From the normalization at infinity, we get that det(TV) = 1. □ 

8.2. Local Analysis Near Univalent Ends. In the previous section we described the main term 
of the asymptotics of S away from A. In this section we shall do the same near the points in A. 
Recall that there exists exactly one k — k(a) such that the arc A^ is incident with a. Until the end 
of this section, we understand that k is this fixed integer. Moreover, we let J a to be the possibly 
empty index set such that A° has a as its endpoint for each j £ J a . 

8.2.1. Riemann- Hilbert Problem for Local Parametrix. Let 5 > be small enough so that the 
intersection of the ball of radius S centered at a, say U a ,s, with each of the arcs comprising £ and 
incident with a is again a Jordan arc. We are seeking the solution of the following RHP-Pa : 

(a) V a is analytic in U a ,s \ S; 

(b) V a has continuous traces on each side of U a .s H £ that satisfy RHP-6>(bl)-(b3); 

(c) V a has the behavior near a within U a .s described by RHP-6>(c); 

(d) VaAf- 1 =1 + 0(l/n) uniformly on dU a ,s \ S, where TV is given by (8.1). 
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We solve RHP-'P a only for n E N e . For these indices the above problem is well-posed as det(7V) = 1 
by Lemma 4 and therefore A/" -1 is an analytic matrix function in D a . In fact, the solution does 
not depend on the actual value of e, however, the term C(l/n) in RHP-'Pa(d) does depend on e as 
well as 5. That is, this estimate is uniform with n for each fixed e and 5, but is not uniform with 
respect to e or S approaching zero. 

To describe the solution of RHP-Pa, we need to define three special objects. The first one is 
the so-called G-function whose square conformally maps U ai s into some neighborhood of the origin 
in such a fashion that Afc is mapped into negative reals. The second one is a holomorphic matrix 
function needed to satisfy RHP-'Pa (d). The third is a holomorphic matrix function that solves 
auxiliary Riemann-Hilbcrt problem with constant jumps. 

8.2.2. G -Function. Set 

G a (z) := f h(t)dt, zeU aJ \A k . 

J a 

Then G a is a holomorphic function in U a ,s \ Afc such that 

(8.3) |$e- G »|=l in U a , s - 

Indeed, since both a\ and a belong to A and the Green differential dG has purely imaginary 
periods, the integral J dG is purely imaginary itself. It is also true that 

(8.4) G++G- = on A fe n U a , s 

since h + + h~ = on A. Moreover, it holds that the traces G^ have purely imaginary values on 
Afc as the same is true for h ± (t)dt (recall that the quadratic differential h 2 (z)dz 2 is negative on 
Afc). The last observation and (8.4) imply that G\ is a holomorphic function in U a .s that assumes 
negative values on Afc n U a ,s- Furthermore, 

(8.5) |(G^(a)|=21im|/ l (t)||t-a| 1 /V0. 

Property (8.5) implies that G\ is univalent in some neighborhood of a. Without loss of generality, 
we can assume that 8 is small enough for G\ to be univalent in U a ,s- Hence, G\ maps U a ,& 
conformally onto some neighborhood of the origin. In particular, this means that Afc can be 
extended as an analytic arc beyond a by the preimage of [0, oo) under G\ and we denote by Afc 
this extension. 

Let I + := {z : Arg(z) = 2tt/3}, I := {z : Arg(z) = rr}, and Z_ := {z : Arg(z) = -2rr/3} be 
three semi-infinite rays oriented towards the origin. Since we had some freedom in choosing the 
arcs Afc± , we require that 

G 2 a ((A fc+ U Afc_) n U a ,s) C /+ U /_. 
The latter is possible as G\ is conformal around a. We denote by U^ s (resp. U~ s ) the open 
subset of U a _s that is mapped by G\ into the upper (resp. lower) half-plane. Clearly, there are 
two possibilities, either Afc + C U^ s and therefore Afc is oriented towards a, or Afc + C U~ s and 
respectively Afc is oriented away from a (see Figure 6). 

Finally, since the traces G^ are purely imaginary on Afc n U a j, satisfy (8.4) there, and the 
increment of the argument of G a is 7r when a is encircled in the clockwise direction from a point 
on Afc n U a ,s back to itself, we can define the square root of G a that satisfies 

(8.6) (G a /2 ) + = ±i{G)/ 2 y on AfcHC/a,*, 

where the sign + must be used when Afc is oriented towards a and the sign — otherwise. 
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FIGURE 6. Schematic representation of the arcs Afc, A*,, Ak±, the domains U^ s (shaded 
part of the disk) and U~ s (unshaded part of the disk), and two possible cases: Ak± C U^ s 
(Afc oriented towards a) and Afc± C U^ s (Afc oriented away from o). 



8.2.3. Matrix Function £ . Let arg(z — a) be the branch of the argument of (z — a) that was used 
in the definition of p in (2.11). Without loss of generality we assume that its branch cut is Afc. 
Put 

(8.7) W:= ($ n wl /2 (a--) aa/2 y 3 in U a j\A kl 

where we take the principal value of the square root of Wk (we assume that S is small enough so Wk 
is holomorphic and non-vanishing in U a ,s) and use the branch (5.14) to define (a — •) Q °/ 2 . Then it 
holds that 

(8.8) (wl /2 (a-r« /2 ) 2 =e ±a ^p in U± s 
and 

(8.9) (wl /2 ( a --) a ^Y {wl/ 2 (a--) a ^y =p on A fc n U a , s . 
So the matrix function AfW is holomorphic in U a j \ (Afc U [jj A") and 

Ar+W+=M-W-(° on A k nU a ,s- 

Moreover, it is, in fact, holomorphic across each A", j € J a , as 

by (7.4) and since W is diagonal. Hence, we deduce from (8.4) that AfW exp{— nG a a^\ is holo- 
morphic in U a ,s \ Afc and has the same jump across Afc as AfW. Define 

£ :=MWcxp{-nG a a 3 }-^ ( ^ ) (7rnG a ) CT3/2 

where the sign — must be used when Afc is oriented towards a and the sign + otherwise. Since the 
product 

.(«^)T(^^)(-!i)(4l)((^)T 

is equal to X by (8.6), the matrix function £ is holomorphic in U a< s \ {a}. Now, the second part 
of (8.2) and (8.7) yield that all the entries of AfW behave like 0(\z - a|~ 1/4 ) asz^a. Hence, 
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it follows from (8.5) that the entries of £ can have at most square-root singularity there, which is 
possible only if £ is analytic in the whole disk U a $■ 

8.2.4. Matrix Functions iff and . The following construction was introduced in [32, Theorem 6.3]. 
Let I a and K a be the modified Bessel functions and H<p and be the Hankel functions [1, 
Ch. 9]. Set iff to be the following sectionally holomorphic matrix function: 

/ /42C 1 / 2 ) ^ Q (2C 1/2 ) \ 

~ I 2 m C /2 r a (2C 1/2 ) -2C 1/2 i^ (2C 1 / 2 ) ) 




(2(-C) 1/2 ) \H ( a 2) (2(-C) 1 /2) 

e' 2 (2(-C) 1/2 ) K 1/2 (hL 2) )' (2(~C) 1/2 



\H {2) (2(-C) 1 / 2 ) -| J ffi 1) (2(-C) 1 / 2 ) 

-ttC 1 / 2 (i/i 2) )' (2(-C) 1/2 ) K 1/2 (l^)' (2(-C) 1/2 ) 

for — 7r < Arg(C) < — 27r/3, where Arg(£) G (— 7T, 7t] is the principal determination of the argument 
of Q. Assume that the rays /, 7+, and /_ defined in Section 8.2.2 are oriented towards the origin. 
Using known properties of I a , K a , Ha \ Ha , and their derivatives, it can be checked that ^ is 
the solution of the following Riemann-Hilbert problem RHP-^: 

(a) \t is a holomorphic matrix function in C \ (/ U i+ U i_); 

(b) iff has continuous traces on I + U I_ U / that satisfy 



O 
O 

o 
o 



if> + = 



(c) iff has the following behavior near 0: 

| C |a/2 | C | Q /2 
| C |a/2 | C |a/2 

log ICI log ICI 
log ICI ICI 

id Q/2 icr a/2 
ici q/2 icr a/2 

| C |- Q /2 | C |-a/a 
| C |-a/a | C |-a/2 

(d) iff has the following behavior near 00 



1 

1 

-1 







on I± 
on /: 



if a < 0, as ( — » 0, 

if a = 0, as ( — » 0, 

if a>0, as C^0 in |Arg(C)| < 2tt/3, 

if a>0, as ( -> in 2tt/3 < |Arg(C)| < tt. 



*(C) = (2K 1 / 2 ) 



-V3/2 1 



1 i 
i 1 



I + 0(r 1 / 2 ))cxp{2C 1 / 2 ( T3} 



uniformly in C \ (JU 1+ U J_). 

Finally, if we set iff := a 3 ^>o- 3 . It can be readily checked that this matrix function satisfies 
RHP-^ with the orientations of the rays /, I + , and /_ reversed. 
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8.2.5. Solution of RHP-TV With the notation introduced above, the following lemma holds. 
Lemma 5. For n £ N E , a solution of RHP-'P Q is given by 

(8.10) V a =£^W-\ (=(n/2) 2 G 2 a , 

if Afe is oriented towards a and with ^ replaced by \& otherwise, where ^ = $(-;a a ). 

Proof. Assume that A^, and respectively A k ±, is oriented towards a. In this case G 2 preserves the 
orientation of these arcs and we use (8.10) with The analyticity of £ implies that the jumps of 
V a are those of ^W -1 . By the very definition of G 2 and the latter has jumps only on S n U a ,s 
and otherwise is holomorphic. This shows the validity of RHP-TV a )- It also can be readily verified 
that RHP-Pa (b) is fulfilled by using (7.4), (8.8), and (8.9). Next, observe that RHP-7V(c) follows 
from RHP-I'(c) upon recalling that \G 2 a (z)\ ~ \z - a\ and \W{z)\ ~ \z - a | (Q » /2)<T3 as z -> a. 
Observe now that with ( defined as in (8.10), it holds by the definition of £ and RHP-^E'(d) that 

TaAf- 1 -1 = JVW exp {-nG a cr 3 } O (~\ exp {nG a a 3 } W^AV 1 = MO f ~) M' 1 

on dU at s, where we also used (8.3). Multiplying the last three matrices out we get that the entires 
of thus obtained matrix contain all possible products of S n / S n (oo), hS*/ S n (oo), S n -i/S*_i(oo), 
and /i«5*_ 1 /5*_ 1 (oo). Then it follows from (6.11) used with e < 5 that the moduli of the entires of 
V a M~ x — T are of order O (1/n) uniformly for n € N E . This finishes the proof of the lemma since 
the case where is oriented away from a can be examined analogously. □ 

8.3. Local Analysis Near Trivalent Ends. In this section we continue to investigate the be- 
havior of S near the points in E. However, now we concentrate on the zeros of h, that is, the 
trivalent ends of A. As in the construction of the lens E, let Af,^ be the arcs comprising A incident 
with b which are numbered in the counter-clockwise fashion. 

8.3.1. Riemann- Hilbert Problem for Local Parametrix. As before, we denote by Ubfi a disk centered 
at b of small enough radius 5 ( "small enough" is specified as we proceed with the solution of RHP- 
Vh). We are seeking the solution of the following RHP-TV 

(a) Vb is analytic in Ub.& \ E; 

(b) Vb has continuous traces on each side of Ub,s H E that satisfy RHP-5(b2) and (b4); 

(c) Vb is bounded in the vicinity of b; 

(d) VbAf- 1 =1 + 0(l/n) uniformly on dU b j \ E, where J\f is given by (8.1). 

As in the case of RHP-TV we consider only the indices n £ N £ and the estimate 0(1/ n) in 
RHP-"Pf,(d) is not uniform with respect to s or 8 approaching zero. 

8.3.2. G-Function. Set, for convenience, Sb,k to be the sectorial subset of Ub,s bounded by A^fc+i, 
A^fc-i, and dUb : s- Define 

G b (z) := {-l) k f h(t)dt, z £ S btk . 

Jb 

Thus defined, the function Gb satisfies 

(8.11) |$6 Gb | = l in S b ,iUS b ,3 and |$^ Gb |=l in S b>2 - 

The same reasoning as in (8.4) yields that Gb is a holomorphic function in Ub.s \ Af,^ whose traces 
on Abj n Ub,s and Afc i3 n Ub,s as well as on both side of A b ,2 H U b ,s are purely imaginary and 
G^ + G^ = on Ab^ H tig. The last observation implies that G\ is a holomorphic function in Ub : $ 
for S small enough that assumes negative values on each A^ n Ub : s- 
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Recall that h 2 has a simple zero at b and therefore \h(z)\ 



&I 1 / 2 as z — > b. This, in turn, 



implies we can holomorphically define a cubic root of G\ in Ub,s- In what follows, we set G 



2/3 



to 



be a conformal map of Ub y s onto some neighborhood of the origin such that 
(8.12) G^ /3 (A b , fc n[/ M )c/ fe , 

where we set Ik := {z : arg(z) = ir(2k — l)/3}, k G {1,2,3}, and these rays are oriented towards 
the origin. Moreover, since we had some freedom in choosing the arcs A^, we require that 



(8.13) 



G 



2/3 



A 6ife nc/ M c h, 



where Ik := {z : Arg(z) = 27r(fc — 2)/3}, k <E {1, 2, 3}, and the rays are once again oriented towards 
the origin. Such a choice is always possible as G\ maps the sector Sb t k 3 Af,^ onto a neighborhood 
of the origin cut along 1% . 



Finally, since the traces G b are purely imaginary on A^ 2 H Ub,s, satisfy G b 



-G b there, and 



the increment of the argument of Gb is 3n when b is encircled in the clockwise direction from a 
point on A&.2 H Ub,s back to itself, we can define the sixth root of Gb, which is equivalent to the 

2/3 

fourth root of GJ , that satisfies 



(8.14) 



(gT) + =±*(gT) 



on A ij2 n 



where the + sign must be used when Af,^ are oriented towards b and the — sign otherwise. 

8.3.3. Matrix Function £. Recall our notation, pb.k = P\A b k - For each function pb,k, fix a contin- 
uous branch of the square root yjpb,k- Let TV be the solution of RHP-TV presented in Section 8.1. 
Set 



W := < 



\fPbJ 



in S b ,i, 



in Sb.2 U Si 



'6,3- 



Then the matrix function TVW is holomorphic in Ub,s \ UfcA^fc and satisfies 



TV+W+ =TV_W_ 













on 








K-J 


») 


on 



on (A 6 , 2 U A M ) n £/&„ 



Further, put 



TVyVexp{nG 6 cr 3 } . 



£* := 



TVW 



in SVi U 5^3, 
exp{7iG h cr 3 } , in S b . 2 , 



±1 
V Tl 

where we use the upper signs when the arcs A& are oriented towards b and the lower ones otherwise. 
Then £* is a holomorphic matrix function in Ub,s \ A^ and 

1 



since Gt 



£% = £*_ 



-G b there. Finally, define 



-1 



£ :=£* 



1/2 T l/2 
Ti/2 -i/2 



on A b .2 n U b ,s 



G£ 3/6 (3n/2) CT3 / 6 
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where we use the — sign when the arcs are oriented towards b and the + sign otherwise. Since 
the product 



{Gl /e ) 



1 ±i 



1 
-1 



1/2 T l/2 
=F*/2 -i/2 



is equal to I by (8.14) where we use the upper signs when the arcs are oriented towards b and the 
lower ones otherwise, the matrix function £ is holomorphic in Ub,s\{b}. Since |Gb(^)| 1 / 6 ~ \z — b] 1 / 4 
as z — > b and by the first part of (8.2), the entries of £ can have at most square-root singularity 
there. Therefore £ is analytic in the whole disk Uf, g. 

8.3.4. Matrix Function T. The following construction is a modification of the one introduced in 
[15, Section 7]. Let Ai be the Airy function. Set 



and 

Further, put 
T := < 



T (C) 



Ti(C) 



Ai(C) Ai(e^c) 
Ai'(C) e^Ai'fe^C 




-e 4 3* Ai [e 2 " £ 
-Ai' (e^C 



e t 



e e 



To, 
To 



1 
-1 1 

It is known that 



Arg(C) G (0, f ) 
Arg(C)e (f,f 
Arg(C)e (f ,tt 



■= < T 1; 
Ti 



- 
1 

1 
1 1 



Arg(C)e (-f,0), 
Arg(C)e (-f ,-f) 



Arg(C) G (-7T, - 



271 \ 



Ai(o = ^r 1/4 ex P {-ic 3 / 2 }(i + o(r 3/2 )) 
Ai '(0 = -2^c 1/4 ex P {-ic 3 / 2 }(i + o(r 3/2 )) 

as C — » oo in the angle |Arg(£)| < 7r, from which it was deduced in [15, Lemma 7.4] that 



(8.15) 



T(C) 



2^ 



-cr 3 /4 



1 i 

-1 i 



Z + O 



exp 



-3C ^ 



as C -> 00 for |Arg(C)| G (§, x) U tt) . Asymptotics for |Arg(C)| G (0, §) can be obtained 
by multiplying the left-hand side of (8.15) by the matrix 



1 







from the right. Altogether, 



it can be checked as in [15, Section 7] that T is the solution of the following Ricmann-Hilbcrt 
problem RHP-T: 

(a) T is a holomorphic matrix function in C \ {J k (jk U I^j ; 

(b) T has continuous traces on (J fc I Ik U Ik ) that satisfy the jump relations described by 



Figure 7; 

(c) each entry of T has a finite nonzero limit at the origin from within each sector; 

(d) the behavior of T near 00 is governed by (8.15). 
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FIGURE 7. The jump matrices that describe the relations between the traces of T on 



Finally, if we set T := (73X0-3. It can be readily checked that this matrix function satisfies 
RHP-T with the orientations of the rays Ik and Ik reversed. 

8.3.5. Solution of KHP-Vb- With the notation introduced above, the following lemma holds. 

Lemma 6. For n € N e , a solution of RHP-T 5 ;, is given by 

(8.16) V b = 2V^eT£TW-\ C = (3n/2) 2 / 3 G^ /3 , 

if the arcs A^j. are oriented towards b and with T replaced by T otherwise. 

Proof. Assume that the arcs Ab,k arc oriented towards b. As £ is holomorphic in Ub,s, it can be 
readily verified using (8.12) and (8.13) that TW _1 satisfies RHP-T^b). It is also evident that 
TW _1 has no other jumps and therefore RHP-T^a) is fulfilled. Since all the matrices are bounded 
in the vicinity of 6, so is RHP-'Pb(c). Observe now that with Q defined as in (8.16), it holds by the 
definition of £ and (8.15) that 

VbAf^ 1 -1 = JVW cxp {nG b o- 3 } ° (j^j CX P {-nGWs} WW -1 

on dU b ,8 n d(S b ,i U S b ,a) and 

VbAf- 1 -I = jm(° 1 J J cxp {nG b a 3 } O cxp {-nG b <J 3 } ( ° ~J J W^AT 1 

on dUb,s r\dSb : 2- As in the case of RHP-'P a (d), these representations yield RHP-T^d) on account 
of (8.11) and (6.11) used with e < 6. This finishes the proof of the lemma since the case where the 
arcs Aft.fc are oriented away from b can be examined analogously. □ 

8.4. Final Transformation. Denote by E the reduced system of contours that we define as 

9 



E:= E\ 



U (J dU e , s 



A U |J U e , s U (J A.% 

V L e£E fc=l 

(see Figure 8). For this new system we consider the following Riemann-Hilbert problem (RHP-7?.): 
(a) 72. is a holomorphic matrix function in C \ E and TZ(oo) = I; 
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A 5 



A 4 




Ai 



FIGURE 8. Contour E for £ from Figure 5. 



(b) the traces of 1Z on each side of E are continuous except for the branching points of E, where 
they have definite limits from each sector and along each Jordan arc in E. Moreover, they 
satisfy 



uniformly on IJeeE dU e ^- Further, as $ _2n converges to zero geometrically fast away from A, the 
jump of 1Z on E \ [J eeE dU e ,s is geometrically uniformly close to I. Hence, (8.18) holds uniformly 
on E. Thus, by [14, Corollary 7.108], RHP-7?. is solvable for all n large enough and 1Z± converge to 
zero on E in L 2 -sense as fast as 1/n. The latter yields (8.17) locally uniformly in C \ E. To show 
that (8.17) holds at z g E, deform E to a new contour that avoids z. As the jump in TLHP-1Z is 
given by analytic matrix functions, one can state an equivalent problem on this new contour, the 
solution to which is an analytic continuation of 71. However, now we have that (8.17) holds locally 
around z. Compactness of E finishes the proof of (8.17). □ 

Now, it can be verified directly from Lemmas 4, 5, and 7 that the following lemma holds. 

Lemma 8. The solution of RHP-5 exists for all n S N e large enough and is given by 



\ lZV e , in U e .s for each ee£, 
where 1Z is the solution of RHP -7Z. 

It is an immediate consequence of Lemmas 2, 3, and 8 that the following result holds. 




on dU C: s for each e G E, 




Then the following lemma takes place. 

Lemma 7. The solution of RHP -TZ exists for all n£N E large enough and satisfies 

(8.17) K=l + 0(l/n), 

where 0(\/n) holds uniformly in C (but not uniformly with respect to e). 

Proof. By RHP--p a (d) and RHP--p b (d), we have that RHP-ft(b) can be written as 

(8.18) 1Z+ =1l-(X + 0(l/n)) 



(8.19) 
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Lemma 9. If Condition GP is fulfilled, then the solution of RHP-} 7 uniquely exists for all n € N e 
large enough and can be expressed by reversing the transformations y — > T — > S using (7.2) and 
(7.6) with S given by (8.19). 

9. ASYMPTOTICS OF NUTTALL-STAHL POLYNOMIALS 

9.1. Proof of Theorem 4. Assume that n £ N £ . For any given closed set in D* , it can be easily 
arranged that this set lies exterior to the lens S. Thus, the matrix y on this closed set is given by 

where TZ is the solution of RHP-1Z given by Lemma 8 and N is the solution of RHP-A/" given by 
(8.1). Then 



(9.1) KAf = 



(1 + fnl) g 7— s +^n2-j-^ 7—T {l+V n l)-=-T— r+V n2 — -— 

SW(oo) QS^lpo) S„{oo) S£_i(oo) 

\ u «3^-7 — r + (1 + v ni ) — — r w„ 3 -?T7 — \ + \ l + v ni — 7 — r- / 



with |w„fe| < c(e)/n uniformly in C by (8.17); and therefore (2.28) follow from (7.1). 

To derive asymptotic behavior of q n and R n on A\£, we need to consider what happens within 
the lens £ and outside the disks U e .s- We shall consider the asymptotics of y from within [j. ilk+, 
the "upper" part of the lens £, the behavior of y in (J fc Qk- can be derived in a similar fashion. 
We deduce from Lemma 9 that 

y+ = 7 r (*ao + ( (4+) L /p ; ) (* + r- = ir (^o + ( ( ($+ % 

locally uniformly on A\E. Therefore, it holds that 

= (^AA)+ 1 (7A$ + )" + (^)+72( < f + )"7P, 

y+ = (7^)+( 7 a/$ + )". 
Hence, we get (2.29) from (9.1), (7.1), and (2.25). □ 

9.2. Proof of Corollary 5. Since p — [n/n]p = R n /q n , it follows from (2.28) that 

~_r , , = S* h 1 + Vnl + K 2 ^)(7^/7n)(S*-l/S*) 

P i n l n \ P Sn $2 „ 1 + Vni + ( Un2 /$)( 7 */ 7T[ )(5 n _ l/l5n ) ■ 

Since v n i is vanishing at infinity, it follows from (8.17) that |u„2 < f|, |uri2/ < £ > | — c(s)/n uniformly in 
D*. Thus, (2.30) is the consequence of (2.27). □ 
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